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Abstract 

We review the current theory of how galaxies form within the cosmological framework provided by the cold dark 
matter paradigm for structure formation. Beginning with the pre-galactic evolution of baryonic material we describe 
the analytical and numerical understanding of how baryons condense into galaxies, what determines the structure of 
those galaxies and how internal and external processes (including star formation, merging, active galactic nuclei, etc.) 
determine their gross properties and evolution. Throughout, we highlight successes and failures of current galaxy 
formation theory. We include a review of computational implementations of galaxy formation theory and assess their 
ability to provide reliable modeling of this complex phenomenon. We finish with a discussion of several "hot topics" 
in contemporary galaxy formation theory and assess future directions for this field. 
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1. Introduction 

While the concept of galaxies as "island universes" can be traced back to IWrighll(ll750h and Kant ('1755') the study 
of the formation of galaxies did not begin until after their extra-Galactic status was confirmed by Hubble C19291) . In 
fact, much of the early work on galaxy evolution and formation was driven by the necessity of understanding galaxies 
in order to answer questions of cosmology (such as whether or not the Universe began with a Big Bang). While an 
understanding of galaxies remains necessary for such reasons even today, the field has since become an important one 
in its own right. 

Modern galaxy formation theory therefore grew out of early studies of cosmology and structure formation and 
is set within the cold dark matter cosmological model and so proceeds via a fundamentally hierarchical paradigm. 
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Observational evidence and theoretical expectations indicate that galaxy formation is an ongoing process which has 
been occurring over the vast majority of the Universe's history. The goal of galaxy formation theory then is to describe 
how simple physics gives rise to the complicated set of phenomena which galaxies encompass. 

Galaxy formation is very much an observationally driven field in the sense that we are still decidedly in the stage 
of making new experimental discoveries rather than performing precision tests of well-specified theoretical models. 
While this situation shows signs of a gradual shift to the "precision tests" phase it seems unlikely that the transition 
will be completed any time soon. In addition, astronomy is perhaps uniquely hindered by experimental biases, since 
we are not able to design the experiment, merely to observe what the Universe has decided to put on show. The 
complicated nature of the resulting selection effects result in a secondary, but very important, role for theoretical 
models, namely in quantifying these biases and interpreting the data. While this secondary role is well established 
it needs to become more so, in particular it should become an integral part of any observational campaign and will 
require direct and simple access to modeling capabilities for all. 



2. Background Material 



Modern galaxy formation theory is set within the larger scale cold dark matter cosmological model ([Blumenthal et al, 



19841) . The success of t hat model in explaining the cosm i c microwave background (CMB') (IKomatsu et i 
and large scale structure ( Seliak et al. , 2005 : Percival et al. , 2007a : Ferramacho et al. , 2008 : Sanchez et all 



1.1 l2009h 



2009) of 



the Universe makes it the de facto standard. However, it is important to recognize that the scales on which the 
simplest cold dark matter (CDM) model has been most precisely tested are much larger than those that matter for 
galaxy formatiorQ. As such, we cannot fully rule out that dark matter is warm or self-interacting, although good 
constraints exist on both of these properties ([M arkevitch et al., 2004; Boehm and Schaeffer, 2005; Ahn and Sh apiroi 
20051 iMiranda and MacciSibOOTtkandall etal.L 120081: 1 Yukselet al.L 120081: iBovarskv et al.ll2008l) . 



Mo re extensive re vie ws of the cold dark matter cosmological model can be found in, for example. lNarlikar and Padmanabhan 



( l2001h . lFrenkl ( 120021) and lBertone et al.. (.2005.) . 



2.1. Background Cosmology 

A combination of experimental measure s , including studies of the CMB dDunklev et al. I l2009h . large scale struc- 
ture ( Tegmark et al. 2004; Cole et al. 20051: Tegrnark et al. 20061 : Percival et al. 20073 ^). t he Type la supernovae 
magnitude- redshift relation ("Kowalski et al."2008') and galaxy clusters ("Mant z et airEo09l: IVikhUnin et al.[ . 120091: 



Rozo et al 



h ave now placed stron g constraints on the parameters of the cold dark matter cosmogony. The 



picture that emerges ( Komatsu et al. , 2010l) is one in which the energy density of the Universe is shared between dark 



energy (Qa = 0.728+° °|^), dark matter (Qc = 0.227 + 0.014) and baryonic matter (Qb = 0.0456 + 0.0016), with a 
Hubble parameter of 70.4;^[^ km/s/Mpc. Perturbations on the uniform model seem to be well described by a scale-free 
primordial power spectrum with power-law index Hs - 0.963 ± 0.012 and amplitude erg = 0.809 ± 0.024. 

Given such a cosmological model, the Universe is 13.75 ± 0.11 Gyr old. Galaxies probably began forming at 
z ~ 20 - 50 w hen the first sufficiently deep dark ma tter potential wells formed to allow gas to cool and condense to 
form galaxies ( Tegmark et al Jl 1 997 : Gao et al ] |2007l) . 



2.2. Structure Formation 

The formation of structure in the Universe is seeded by minute perturbations in matter density expanded to cosmo- 
logical scales by inflation. The dark matter component, having no pressure, must undergo gravitational collapse and, 
as such, these perturbations will grow. The linear theory of cosmological perturbations is well understood and pro- 
vides an accurate description of the early evolution of these perturbations. Once the perturbations become nonlinear, 
their evolution is significantly more complicated, b ut simple arguments (e.g. spherical top-hat collapse and devel- 
opments thereof; Gunn 1977 : Shaw and Motal 2008i) provide insight into the basic behavior. Empirical methods to 



' Although measurements of the Lyman-o; forest iSlosar et~al l l2007l : IViel et alll2008ll and weak-lensing iMandelbaum et alil200d) give inter- 
esting constraints on the distribution of dark matter on small scales they typically require modeling of either the nonlinear evolution of dark matter 
or the behavior of baryons (or both) which complicate any interpretation. They are therefore not as "clean" as CMB and large scale structure 
constraints. 
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determine the statistical distribution of matter i n the nonlinear regime exist ( [Hamilton et all 1991 ; Peacock and Dodds 



1996; Smith et al. 2003; Heitmann et al. 2009). These, together with N-body simulations (e.g. iKlvpin and Shandarin 



1983; Springel et al.. 2005b; Heitmann et al. 200&) show that a network of halos strung along walls and filaments 
forms, creating a cosmic web. This web is consistent with measurements of galaxy and quasar clustering on a wide 
range of scales. 



2.3. Halo Formation 

The final result of the nonlinear evolution of a dark matter density perturbation is the formation of a dark matter 
halo: an approximately stable, near-equilibrium state supported against its own self-gravity by the random motions of 
its constituent particles. In a hierarchical universe the first halos to form will do so from fluctuations on the smallest 
scales. Later generations of halos can be thought of as forming from the merging of these earlier generations of halos. 
For the purposes of galaxy formation, two fundamental properties of the dark matter halos are of primary concern: 
(i) the distribution of their masses at any given redshift and (ii) the distribution of their formation histories (i.e. the 
statistical properties of the halos from which they formed). 



2.3.1. Halo Mass Distribution 



The insight of iPress and Schechten (119741) was that halos could be associated with peaks in the Gaussian random 
density field of dark matter in the early universe. Using the relatively simple statistics of Gaussian random fields they 
were able to derive the following form for the distribution of dark matter halo masses such that the number of halos 
per unit volume in the mass range M to M + SM is 6M(dn/dM) wher^ 



dn 

dM ' 



1/2 



X exp 



PO ^c(0 
M2 cr(M) 

6l(t) 



IcrHM) 



diner 



dlnM 



(1) 



where po is the mean density of the Universe, cr(M) is the fractional root variance in the density field smoothed using 
a top-hat filter that conta ins, on average, a mass M, and 5JJ:) is the critical overdensity for spherical top-hat collapse 
at time / dEke et al.lll996l). 



While the lPress an d Schechten (119741) expression is remarkably accurate given its simplicity, it does not provide a 
sufficiently accura te description of modern N -body measures of the halo mass function. Several attempt s have been 
made to "fix" the IPress and Schechten (119741) theory by appeahng to different filters and barriers (e.g. ISheth et al 



20011) although to date none are able to accurately predict the measured form without adjusting tunable parameters. 
The most accurate fit ting formulae currently available are those of ,Tinker et al.. (,2008; see also ,Reed et al.i i2007t 
Robertson et al ] l2009l) . Specifically, the mass function is given by 



dn _ p d In cr ' 
dM ~ ^ M~dM~' 



(2) 



where 



f(cr) = A 



(3) 



and where A, a, b and c are parameters determined by fitting to the results of N-body simulations. The mass variance 
cr^(M) is determined from the power spectrum of density fluctuations 



Cr2(M) 



^ j| P{k)T\k)Wl{k)k'dk 



(4) 



^The original derivation bv lPress and Schechtej il974) differed by a factor of 2, resulting in only half of the mass of the Universe being locked 
up in hal os. Later de rivations placed th e method on a firmer mathematical basis and resolved this problem, a symptom of the "cloud-in-cloud" 
problem teond et alJl991. : .Bowei„199ldLacev and Coldll993l) . 
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where P(k) is the primordial power sp ectrum (usually taken to b e a scale-free power spectrum with index ris), T(k) is 
the cold dark matter transf er function ( Eisenst ein and Hull 199 9*) and WM{k) is the Fourier transform of the real-space 
top-hat window function. iTinker et al. ( 2008) give values for these parameters as a function of the overdensity. A, 
used to define what we consider to be a halo. Additionally, they find that the best fit values are functions of redshift. 
Earlier studies had hoped to find a "universal" form for the mass function (such that the functional form was always 
the same when expressed in terms of v = cr/6c)- While this is approximately true, the work of .Tinker et al.l (l2008l) 
demonstrates that universality does not hold when high precision results are considered. 



2.3.2. Halo Formation Distribution 

A statistical description of the formation of halos, specifically the sequence of mergi ng events and the masses of 
halos inv olved in those events, c an be extracted using similar arguments as the original Press and Schechtei ( 1974 ) 
approach dLacev and Colelll993l ). These show that the distribution of halo progenitor masses, Mjm at redshift zi for a 
halo of mass M2 at later redshift Z2 is given by 



dA^ 



1/2 



dlncr 
dlnMi 



M2 



X exp 



(5cl - 5c2)' 



where a-\ - (t(M\), 0-2 - cr(M2), dd = <5c(zi), 5c2 = Sdzi)- N Vith a zero time-lag (i.e. 
del — > 6c2) this can be interpreted as a merging rate (although see^ Neistein and Dekell ( 200 . 

Repeated application of this merging rate c an be used to build a merger tree. Finding a suitable algorithm is non 



(5) 

as zi Z2 and therefore 
for a counter argument). 



trivia l and many attempts h ave been made ( Kauffmann and White 1993 ; Somerville and Primack 1999t Cole et al . 
I2OOO';' Parki nson et al ] l2008l) . A recent examination of alternative algorithms is given bv' Zhang et al. Current 



implementations of merger t ree algorithms are high ly accurate and can reproduce the progenitor halo mass distribution 

over large spans of redshift ( Parkinson et al.ll2008 ). 

A fundamental limitation of anv lPress and Schechtei (I1974') based approach is that the merger rates are not sym- 
metric, in the sense that switching th e masses Mi and M2 resul ts in two differ ent predictions for the rate of mergers 
bet ween halos of mass M i and M2. Benson et al. ( 2005b and Benson (20081 showed that a symmetrized form of 
the Parkinson et al. (2008) merger rate function could be made to approximately solve Smoluchowski's coagulation 
equation (Smoluchowksi 1916) and thereby provide a solution free from am biguities. Other empirical determinations 
of merger rates have been made (ICole et al ]l2008HFakhouri and Mall2008allbl ). 

In addition to these purely analytic approaches, numerical studies u tilizing N-body simulations have lead to the 



development of an empiric al understanding of halo formation histories ( Wechsler et al. , 20021: 



and h alo-halo merger rates ( Fakhouri and Ma . 2008aL 2009 : Stewart et al. . 2009; Fak houri et al. 



2010h Merger trees can also be extracted directly from N-body simulations (e.g. Hellv et al 



van den Boschl 12002 ) 
.bOlOtlHopkinset all 
2003a; 'Springel et all 



2005bl) which sidesteps these problems but incorporates any limitations of the simulation (spatial, mass and time res- 
olution), and additionally provides information on the spatial distribution of halos. Such N-body merger trees also 



serve to highlight some, perhaps fundamental, limitations of the'P ress and Schechtei ( 19741 ) type approach. For exam- 
ple, halos in N-body simulations undergo periods of mass loss, which is not expected in pure coagulation scenarios. 
The existence of systems of substructures withi n halos (see Sect ion refsec:galaxyOrbits) can even lead to three-body 
encounters which cause subhalos to be ejected (ISales et al ] l2007h . 



2.4. Halo Structure 

Dark matter halos are ch aracterized by th eir large overdensity with respect to the background density. Spherical 
top-hat collapse models (e.g. Eke et al.lll996h show that the overdensity. A, of a just collapsed halo is approximately 
200, with some dependence on cosmological parameters. This overdensity corresponds, approximately, to the virial- 
ized region of a halo, and allows us to define a virial radius as 



3M 
47rp()A 



1/3 



(6) 
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Studies utilizing N-body simulations show that halos approximately obey the virial theorem within this radius and 
that Tv is a characteristic radi us for the transition fro m ordered inflow of matter (on larger scales) to virialized random 
motions (on smaller scales) (ICole and Lacev Il996h . Dark matter halos have triaxia l shapes and the distribution of 
axial ratios has been well characterized from N-body simulations ( Frenk et al. 1988 ; Jing and Sutd 2002 ; Bett et al 
20071) 



Recent N-body studies jNavarro et al ] |2004HMerritt et al.ll2005l: iPrada et al ]l2006l) indicate that density profiles of 
dark m atter halos hav e an approximately universal form in CDM co smologies that is better des cribed by the Einasto 
profile teinastoll 19651) than the Navarro-Frenk- White (NFW) profile (iNavarro et alJll996[[l997l) . The Einasto density 
profile is given by 

p(r)=p_2exp (--[(—) -ll), (7) 



-J -'!)• 

where r_2 is a characteristic radius at which the logarithmic slope of the density profile equals -2 and a is a parameter 
which controls how rapidly the logarithmic slope varies with radius. The parameter a seems to correla te well with the 
height of the peak from which a halo formed, v = 6c(z)/cr(M) as has been shown bv lGao et al.l (120081) who provide a 
fitting formula, 

0.155 + 0.0095v2 ify< 3.907 



a - 



0.3 



ifv> 3.907, 



which is a good match to halos in the Millennium SimulatiorH dSpringel et alJ2005bl) . The value of r_2 for each halo is 
determined from the known virial radius, r^, and the concentration, c_2 = r^lr-i- The NFW profile has a significantly 
simpler form and is good to within 10-20% making it still useful. It is given by 



p(r) = 4 



Ps 



(r/r,)[l+r/rs]2' 



(9) 



where is a characteristic scale radius and is the density at r = r^. For NFW halos, the concentration is defined as 

cnfw - r^lr^- 

Concentrations are found to depend weakly on halo mass and on redshift and can be predicted from the formation 
history of a halo ( Wechsler et al. 2002). Simple algorithms to approximately determine concentrations have been 
proposed bv INavarro et alJ (Il997l) . lEke et al.l (1200 1[) and Bullock et al. (2001b). More accurate power-law fits have 
also been determined from N-body simulations ("G ao et al.<.2008;iZhao et al.ii2009n . 

Integrals over the density and mass distribution are needed to compute the enclosed mass, velocity dispersion, 
gravitational energy and so on for the halo density profile. For NFW halos the integrals are mostly straightforward, 
although some require numerical calculation. For the Einasto profile some of these may be expressed analytically 
in terms of incomplete gamma functions ( Cardo ne et al.ll2005h . Specifically, expressions for the mass and gravita- 
tional potential are provided by Cardone et al.l (l2005h . other integrals (e.g. gravitational energy) must be computed 
numerically. 



3. Pre-Galactic Evolution of Baryons 

Baryons are initiall y distributed near uniformly — they are expec ted to trace the dark matter distribution on scales 
above the Jeans length ( Arons and Silk|[l968 : Gnedin and Huil[l998 ). To form galaxies they must first be concentrated 
by the forces of gravity which are dominated by the distribution of dark matter In particular, we expect that baryons 
will concentrate towards the deep potential wells of dark matter halos. These should therefore be the sites of galaxy 
formation. 



^We have truncated this fit so that a never exceeds 0.3. iGao et all )2008h were not able to probe the behavior of a in the very high v regime. 
Extrapolating their formula to v > 4 is not justified and we instead choose to truncate it at a maximum of o- = 0.3. 
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3.1. Cold or Hot Accretion (Shocks) 

Baryonic material will be dragged along by the gravitationally dominant dark matter such that dark matter halos 
are expected to accrete baryonic material. How much baryonic material they accrete depends upon the depth of their 
potential well an d the pressure of the bar yons. 

According to Okamoto et al. (l2008al) . the mass of baryons which accrete from the intergalactic medium into a 
galaxy halo during some interval 6t after the Universe has been reionized (i.e. the hydrogen content of the Universe 
has been almost fully ionized as a result of emission from stars and active galactic nuclei (AGN)) is given by 



where M^^c is given below, 



Mb = + Ma, 



K = yexp(--^)Mb 



(10) 



(11) 



and where the sum is taken over all progenitors of the current halo, and fevp is the timescale for gas to evaporate from 
the progenitor halo and is given by 

^v/^s(^evp) if < T^evp? 



if Tv > T, 



(12) 



evp 5 



where Tevp i s the temperature below which gas will be heated and evaporated from the halo, Cs is the sound speed in 
the halo gas. lOkamoto et al.l (l2008al) compute Tevp by finding the equilibrium temperature of gas at an overdensity of 
Aevp = 10^. The accreted mass. Mace, is given by 



Ob - 




M.. = ^ i^^h-m; ifr™.>raee, 



(13) 



Here, Taee is the larger of Tgq and the temperature of intergalactic medium gas adiabatically compressed to the density 
of accreting gas where Tgq is the equilibrium tem perature at which radiati ve cooling balances photoheating for gas at 
the density expected at the virial radius, for which lOkamoto et al.l (l2008al) use one third of the halo overdensity. 

An accretion shock is a generic expectation whenever the gas accretes supersonically as it will do if the halo 
virial temperature exceeds the tem perature of the accreting gas (Binney 1977a). Models of acc r etion shocks have 
been presented by several authors ( Bertschingei 1985 : Tozzi and NormanlbOOli Ivoit et al. 2003 : Book and Benson 
boio") with the general conclusion that the shock occurs at a radius comparable to (or perhaps slightly larger than) 
the virial radius when cooling times are long compared to dynamical times. In the other limit of short cooling times, 
it has long been understood that the shock must instead form at much sm aller radii, close to the fo rming galaxy 
(IRees and Ostriken,ll977tlWhite and Frenklll99lh . For example, according to lRees and Ostrikeii(ll977l) : "Unless pre- 
galactic clouds collapse in an exceedingly homogeneous fashion, their kinetic energy of infall will be thermalized by 
shocks before collapse has proceeded by more than a factor ~ 2. What happens next depends on the relative value of 
the cooling and collapse timescales. Masses in the range lO'^-lO'^M© cool so efficiently that they always collapse 
at the free-fall rate, and pro bably quickly fragment in to stars. Larger masses may, however, experience a quasistatic 
contraction phase. . . ". Thus. lRees and Ostrikeil(ll977h clearly understood the difference between the rapid inflow and 
hydrostatic cooling regimes, and correctly identified the transition mass, suggesting that this be identified this with 
the characteristic stellar mass of galaxies. Accretion in these two regimes may be expected to result in very different 
spatial and spectral dis tributions of cooling radiation, leading to the possibility of observationally distinguishing the 



two types of accretion (IFardal et al.Ll2001l) . 



The distinc tion between these t wo reg imes has alwa ys been an integra l part o f analytic models of galaxy formation, 
beginning with lRees and Ostriker ( 1977 ). For example. White and Frenk ( 1991 ) introduced a transition between rapid 
and slow cooling regimes at the point where cooling and virial radii become equal, or, equivalently, the point at 
which cooling and dynamical times at the halo virial radius become equal. In the rapid cooling regime, the accretion 
rate of gas into the central galaxy was then determined by the cosmological infall rate, while in the slow cooling 
regime the accretion rate was determined by the cooling time in the gas. Their Figure 2 illustrates that the rapid 
cooling regime will occur in low mass halos and at hig h redshifts. All subsequent semi-analytic models of galaxy 
formation (e.g. iKauffmann et al. I ll99l ICole et al.llI994l) have adopted this prescription, or some variant of it, and it 
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has also been validated by ID hydrodynamical simulations (e.g lForcada-Miro and Whitelll997h . The validity of this 
prescription has been confirmed by studies which compared its predictions for the condensed masses of galaxies with 
those from smoothed p article hydrodynamics simulations acro ss the boundary of the rapid to slow cooling transition 



tnose rrom srnootned particle nyarodynarmcs simulations acro ss tne bounaary or tne rapia to slow cooling transition 
dBenson et allbOOlbtl^shida et al.Ll200llHellv et alll2003bl) . although it should be noted that the accuracy of these 
comparisons is less than that at which semi-analytic models are now being used. 

Recent work ha s once a gain focused on the formation of accretion shocks. Recen t 3D hydrodynamical simulations 
( Fardal et alj|2001 : see also Keres et al]|2005 : Ocvirk et alj|2008 : Keres et al]|2009h . have suggested that a significant 
fraction of gas in low mass galaxies has never been shock heated (at least within regions adequately resolved by the 
simulations; for example numerical simulations may not adequatel y resolve shocks in the radiative regime due to 
artificial viscosity, numeric al diffusion and other numerical artifacts; lAgertz et al.ll2007 ). Motivated by these results 
Birnboim and Dekel ( 2003") developed an analytic treatment of accretion shock stabilit}5 The accretion shock relies 
on the presence of a stable atmosphere of post-shock gas to support itself. If cooling times in the post-shock gas 
are sufficiently short, this atmosphere cools an d collapses and can no long er support the shock. The shock therefore 
shrinks to smaller radii, where it can be stable. iBirnboim and Dekell(l2003l) find that 



prK{T^)lu^l?ml^ < 0.0126, 



(14) 



where p is gas density, A(r) is the cooling function for gas at temperature T, T\ - l\6)^uyksNa is the post-shock 
temperature and u is the infall velocity, all evaluated for the pre-shock gas at radius r, is required for the shock to 
be stabl^ For cosmological halos this implies that shocks can only form close to the virial radius in halos with 
mass greater than IO^'Mq for primordial gas (or around lO'^M© for gas of Solar metallicity). These values are found 
to depend only weakly on redshift and are in good agreeme nt with th e results of hydrodynamical simulations. It is 
crucial to note that these new criteria are equivalent to that of IWhite and F renk ( 1991) up to factors of order unity. 

As a result, in low mass halos gas tends to accrete into halos "colcf ' — never being shock heated to the virial 
temperature and proceeding to flow along filaments towards the center of the halo where it will eventually shockQ- 
Halos which do support shocks close to the virial radius are expected to contain a quasi-hydrostatic atmosphere of 
hot gas. The structure of this atmosphere is determined by t he entropy that the gas gains at the acc retion shock and 
that may be later modified by radiative cooling and feedback ( Voit et al]|2003 : McCarthy et al.ll2007l) . In practice, the 
transition from rapid to slow cooling regimes is not sharp — halos able to support a shock at their virial radius still 
contain some unshocked gas; because the halos retain a memory of past accretion and because cold filaments may 
penetrate through the hot halo. At high redshifts in particular, the "cold" accretion mode may be active even in halos 
whose accretion of gas is primarily via an accretion shock close to the virial radius. 

The consequences o f rapid vs. slow cooling regimes for the properties of the galaxy forming warrant further study. 
As Croton et al. I (l2006h have stressed, the absence of a more detailed treatment of the rapid cooling regime may not 
be important since, by definition the gas accretion rate in sma ll halos is limited by the growth of the halo rather than 
by the system's cooling time. In contrast, iBrooks et al.l (l2009l) demonstrate in hydrodynamical simulations that in the 
rapid cooling regime accreted gas can reach the galaxy more rapidly, by virtue of the fact that it does not have to cool 
but instead merely has to free-fall to the center of the halo (starting with a velocity comparable to the virial velocity). 
This results in earlier star formation than if all gas were assumed to be initially shock heated to the virial temperature 
close to the virial radius of the halo. It is also clear that the situation needs to be carefully reassessed in the presence 
of eff'ective feedback schemes that prevent excessive star formation, particularly in the high redshift universe. 

While hot atmospheres of gas are clearly present in massive systems such as groups and clusters of galaxies 
(where the hot gas is easily detected by virtue of its X-ray emission), observational evidence for hot atmospheres of 
gas arising from cosmological infall surrounding lower mass systems, such as massive but isolated gala xies is sparse. 



Interesting upper limits have b een placed on the X-ray emission from massive, nearby spiral galaxies (IBenson et al 
2000b : Rasmussen et al. I l2009h . and the hot component of the Milky Way's halo is constrained by the ram pressure 



"^Unpublished work bv lForcada-Miio and Whitj i 19971) . also utilizing a ID hydrodynamics code, reached similar conclusions. 

'The left-hand side of this expression is equivalent, to order of magnitude, to ?sc/'cool where 4c is the sound-crossing time in the halo and /cool 
is the cooling time in the post-shock gas. This provides some physical insight into this condition: if the post-shock gas can cool too quickly sound 
waves cannot communicate across the halo a nd thereby form a hydrostatic atmosphere which can support a shock front). 

Or, more likely, warm jBland-Hawthoml . r2009,V 

'while this picture seems reasonable on theoretical grounds, it as yet has little direct observational support jSteidel et alil2010l) . 
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that it exerts on the Magellanic Stream dMoore and Davislll994t iMastropietro et al. I l2005bh . Additionally, ultraviolet 
line detections of the so-called warm-hot intergalactic medium (intergalactic gas at temperatures of 10^-lO^K) show 
that some of this material must lie close to the Milky Way (,Wang et al. .2005 



Wilhams et al. 2006). What is known 



is that the Milky \y ay's halo contains a significant mass of cold, neutral g as in the form of high velocity clouds 
( Putman et al. 2003 ). This may indicate, as expected on theoretical grounds ( Grain et al. , 2009t) . that the Milky Way 
is in the transition mass range between purely cold and purely hot accretion. 



3.2. Cooling 

Gas which does experience a strong virial shock will have its kinetic infall energy thermalized and therefore be 
heated to of order the virial temperature 

7 CiM.. tiniM 

(15) 



^ _ 2 GMy /imH 

3 Tv ' 



This gas will proceed to form a hydrostatically supported atmosphere obeying the usual hydrostatic equilibrium equa- 
tion: 

dP GM{r) 

-7- = T^p(r\ (16) 

dr 

where P is the gas pressure, p{r) the gas density and M(r) the total (i.e. dark matter plus baryonic) mass within 
radius r. This distribution may be modified, particularly on small scales, b y other con tributors to supporting the halo 
against gravity, such as turbu lence (Frenk et al. 1999), cosmic ray pressure (iGuo and Oh 2 008) and magnetic pressure 
(iGoncpalves and Friacpal Il999h . In any case, the result ng hot atmosphere will fill the dark matter halo with a density 
still several orders of magnitude lower than typical galactic densities. In the absence of any dissipitative process, the 
gas would remain in this state indefinitely. Fortunately, however, gas is able to cool radiatively and so will eventually 
lose energy and, consequently, pressure support, at which point it must fall towards the center of the gravitational 
potential of the dark matter halo thereby increasing its density. Gas which does not experience a shock close to the 
virial radius can fall almost ballistically towards the halo center. It must still lose its infall energy at some point, 
however, shocking closer to the halo center near the forming galaxy. As such, it too will eventually be heated and 
must cool down (although it will presumably do so much more rapidly due to its higher density). 

In the remainder of this subsection we will review the various mechanisms by which such gas cools. 



3.2.1. Atomic 

For metal rich gas, or gas hot enough to begin to coUisionally ionize hydrogen (i.e. T >10"^K), the primary 
cooling mechanisms at low redshifts are a combination of various atomic processes including recombination radiation, 
collisional excitation and subsequent decay and Bremmsstahlung. In the absence of any external radiation field (see 
Section refsec:PhotoHeat), these are all two-body processes and so the cooling rate is expected to scale as the square 
of the density for gas of fixed chemical composition and temperature. It is usual, therefore, to write the cooling rate 
per unit volume of gas as 



i: = n^A(r,Z), 



(17) 



where hh is the number density of hydrogen (both neutral and ionized) and A(r, Z) is the "cooling function" and 
depends on temperature and chemical composition. Typically, the chemical composition is described by a single 
number, the metallicity Z (defined as the mass fraction of elements heavier than helium), and an assumed set of 
abundance ratios (e.g. primordial or Solar). The cooling function can then be found by first solving for the ionization 
state of the gas assuming collisional ionization equilibrium and then summing the cooling rates from the various above 



mentioned cooling mechanisms. Such calculations can be carried out by using, for example. Cloudy (IFerland et al 
Il998i) . Examples of cooling functions computed in this way are shown in Fig.[T] 

In detail, the cooling function depends not just on the overall metallicity, but on the detailed chemical composition 
of the gas. If the abundances of individual elements are known, the corresponding cooling function is easily calcu- 
lated. In theoretical calculations of galaxy formation it is often computationally impractical to follow the abundance of 



numerous chemical species in a large number of galaxies and dark matter halos. Fortunate Iv. lMartmez-Serrano et al 



(120081) have demonstrated that an optimal linear combination of abundances, which minimizes the variance between 
cooling/heating rates computed using that linear combination as a parameter and a full calculation using all abun- 
dances, provides very accurate estimates of cooling rates. The best linear combination turns out to be a function 
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Figure 1 : Cooling functions for gas in collisional ionization equilibrium aie shown as a function of temperature. The solid line corresponds to gas 
of primordial composition, while the dashed line corresponds to gas of Solar composition. 



of temperature. This reduces the problem to tracking the optimal combination of elements for a small number of 
temperatures. This number can then be used in place of metallicity when computing cooling functions. 

Most calculations of cooling rates in cosmological halos assume that the gas is in collisional or photoionization 
equilibrium. Even if the gas begins in such a state, as it cools it can drift away from equilibrium as, particularly at low 
temperatures, the ion-electron recombination timescales can significantly exceed the cooling timescales — as such, the 
ionization state always lags behind the equilibrium state due to the rapidly changing gas temperature. Figure|2]shows 
calculations of effective cooling functions for gas initially in collisional ionization equih brium (and with no externa l 
radiation field) with a fully time-dependent calculation of cooling and ionization state (iGnat and Stern berg '2007V 



Significant differences can be seen, resulting in cooling timescales being a factor of 2-3 longer when non-equilibrium 
effects are taken into account. 

3.2.2. Compton Cooling 

At high redshifts, the density of cosmic microwave background photons becomes sufficiently high that the fre- 
quent Compton scattering of these photons from electrons in the ionized plasma inside dark matter halos results in 
significant cooling of that plasma (assuming that its tempe rature exceeds that of the cosmic microwave background). 
The Compton cooling timescale is given by (|Peebledll968h 



3meC(l + l/Xe) 
TCompton - " -4 " " — — , 

»Crjfl_/pj^g(i - 7CMB/-'ej 

where - n^/nt, rif. is the electron number density, «( is the number density of all atoms and ions, Tcmb is the CMB 
temperature and Tg is the electron temperature of the gas. Unlike the various atomic cooling processes described in 
Section refsec:Cooling Atomic, the Compton cooling rate per unit volume does not scale as the square of the particle 
density, since it involves an interaction between a particle and a cosmic microwave background photon. Since the 
CMB photon density is the same everywhere, independent of the local gas density, the Compton cooling rate scales 
linearly with density. As a result, the Compton cooling timescale is independent of gas density. Thus, if gas in a dark 
matter halo is able to cool via Compton cooling, it can do so at all radii (assuming that it is isothermal and has the 
same electron fraction everywhere). 
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Figure 2: A comparison of cooling functions for gas in collisional ionization equilibrium (upper panel) with effective cooling functions for gas in 
which the time-dependent ionization state is computed throughout the cooling (lower panel). Line types indicate different metallicities (shown in 
the upper panel, values in units of Solar metallicity), while IB and IC labels indicate whether the gas was assumed to be cooling isobarically or 
isochorically. The primary coolant at each temperature is indicated in the upper panel. 
Source: Reproduced with permission from lGnat and Stemberdl2b07l . 
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In a Universe with WMAP 5 yr cosmological parameters, the Compton cooling time is less than the age of the 
Universe above z ~ 6. 



3.2.3. Molecular Hydrogen Cooling 

Somewhat ironically, the gas in dark matter halos at the highest redshifts is too cold to cool any further In halos 
with virial temperatures below lO'^K even shock heated gas (if shocks can occur) will be mostly neutral and therefore 
unable to cool via the usual atomic processes described in Section refsec:CoolingAtomic. Studies show that the 



domi nant coolant in such cases becomes the small fraction of hydrogen in the form of r nolecular hydrog en dAbel et al 
19971). Cooling via molecular hydrogen is crucial for the formation of the first stars dAbel et alj|2002l) and galaxies 



( Bronun et alj|2009l) a nd, therefore, for the sources which cause the reionization of the Universe (iBenson et alj|2006c 
Wise and AbeJ2008al) . The details of molecular hydrogen cooling are more complicated than those of atomic cooling: 
in addition to uncertainties in the molecular chemistry ( Glover and Abell 12008.). i n many cases equilibrium is no t 
reached and the photon background can lead to both negative ( Wise and Abelll2007a ) and positive ( Ricotti et alj2001 ) 

feedbacks. 

Cosmologically, molecular hydrogen forms via the gas-phase reactions ( lMcDowelllll96lh : 



H + e 
H +H 



a ,+y 

H2 +e~ 



(19) 



and 



H-" +H 
+ H 



H2 -hH"^. 



(20) 



Our discussion of molecular hydrogen cooling will mostly follow f^shida et aL I (l2006h . 

The molecular hydrogen cooling timescale is found by first estimating the abundance, /h,,c, of molecular hydrogen 
that would be present if there is no backgroun d of H2-dissociating ra diation from stars. For gas with hydrogen number 
density «h and temperature Ty the fraction is ( Tegmark et alJll997l) : 



3.5 X 10"'*r]-^2[l + (7.4 X 10*^(1 + z?-" 
X exp {-3173/(1 +z))/nHi)]"', 



(21) 



where 7^3 is the temperature Ty in units of lOOOK and «hi is the hydrogen density in units of cm Using this initial 
abundance we calculate the final H2 abundance, still in the absence of a photodissociating background, as 



/H.-/H..exp(^^) 



where the exponential cut-off is included to account for coUisional dissociation of Ho . as inlBenson et al. |(|2006h . 
Finally, the cooling timescale due to molecular hydrogen can be computed using (lGamandPallalll998l) : 



TH, = 6.56419 X 10-^^rv/H>H'iAH;, 



where 



Ah, 



A 



LTE 



1 -1- «cr/nH 

Ah, (LTE) 
Ah,[mh 0] 



(22) 

(23) 

(24) 
(25) 



and 



logio AhJ«h ^ 0] = 



-103 -H 97.59 ln(r) 
-48.05 ln(r)2 -I- 10.8 ln(r)^ 
-0.9032 In(r)'* 



(26) 
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is the cooling functio n in the low density limit (independent of hydrogen density) and we have used the fit given by 
Galh and Pallal ( Il998h . and 

Alte = A,- + A, (27) 
is the cooling function in local thermodynamic equilibrium and 

'9.5 X IQ-^^tP^ 



1 



1 +0.i2r2i 





/ 


0.13 


3\ 


exp 


V 







+3xl0-24exp(- — 

7^3 



ergs cm"' s ' 



A„ 



+ 1.6 X 10"' 



exp - 



11.7 



ergs cm s 



(28) 



(29) 



are the cooling functions for rotational and vibrational transitions in H2 dHoUenbach and McKeelll979l) . 

It is also possible to e stimate the rate of molecular hydrogen formation on dust grains using the approach of 
Cazaux and SpaansI (120041) . who find that the rate of H2 formatio n via th i s rout e can be important in the high redshift 
Universe. In this case we have to modify equation (13) of Tegmark et al.l (Il997h . which gives the rate of change of the 
H2 fraction, to account for the dust grain growth path. The molecular hydrogen fraction growth rate becomes: 

/ = kaf( l-x-2f) + k^n(l-x- 2f)x, (30) 

where / is the fraction of H2 by number, x is the ionization fraction of H which has total number density n. 



ki^ 3.025 X 10-'' -^SniT) 



lOOK 



(31) 



is the dust formation rate coeflicient ( Cazaux and Spaansll2004t eqn. 4), and k ig is the effective rat e coeffi cient for H2 
formation (Tegm ark et al.lll997t eqn. 14). Adopting the expression given bv ICazaux and Spaan s 

(2004 

eqn. 3) for 

the H sticking coefficient, SaiT) and - 0.53Z for the dust-to-gas mass ratio as suggested bv .Cazaux and Spaans 
(120041) . results in * 0.01 for Solar metallicity. This must be solved simultaneously with the recombi nation equation 



govern ing the ionized fraction x. The solution, assuming x{t) - xo/(l +xonkit) and 1 -x-2/ ^ 1 as do lTegmark et al 
jl997l) . is 



fit) = /o-j-exp 



Td 



Tr + t 



(32) 



where Ti- - I /xo/rtu/ki, Td - l/ne/^d, is the hydrogen recombination coefficient and Ei is the exponential integral. 
3.3. Heating 

While gas must cool in order to collapse and form a galaxy, there are several physical processes which instead 
heat the gas. In this subsection we review the nature and effects of those mechanisms. 

3.3.1. Photoheating 

Immediately after the epoch of cosmological recombination, the cosmological background light consists of just 
the blackbody radiation of the CMB. Once stars (and perhaps AGN) begin to form they emit photons over a range of 
energies, including some at energies greater than the ionization edges of important coolants such as hydrogen, helium 
and heavier elements. Such photons can, in principle, photoionize atoms and ions in dark matter halos. This changes 
the ionization balance in the halo and heats the gas (via the excess energy of the photon above the ionization edge), 
thereby altering the rate at which this gas can cool to form a galaxy. 

The ionizing background at some wavelength. A, and redshift, z, is given by 



exp[-T(z,z',/l)] df 



dz' 



(33) 
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Figure 3: Cooling functions for gas at the mean density of cosmological halos at three different redshifts (as indicated in the figure). The dotted line 
shows the cooling function in the absence of any photoionizing backgro und. Ot her lines show the net cooling/heating function in the presence of the 
photoionizing background self-consistently computed by Benson et al.l j20023) . All calculations assume a metallicity of 0.32©. The discontinuity 
at low temperatures shows the transition from net heating to net cooling. 



where £(A,z) is the proper volume-averaged emissivity at redshift z and wavelength A and t(zi,Z2, A) is the optical 
depth between zi and zi for a photon with wavelength A at zi- Calculation of the background requires a knowledge 
of the emissivity history of the Universe and the optical depth (itself a function of the ionization state and density 
distribution in the IGM) as a function of redshift. 

Detailed theoretical models of the ionizing background, using observational constraints on the emissivity history, 
to gether with models of the distributi on of neutral gas and cal culations of radiative transfer, have been developed 
by iHaardt and Madaul (119961: see also iMadau and Haardl I2OO9I) and allow the background to be computed for any 
redshift. 

Once the background is known, its effects on g as in dark matter hal os can be determined. Such calculations require 
the use of photoionization codes, such as Cloudy ( Ferland et al.ll 19981) . to solve the complex set of coupled equations 
that describe the photoionization equilibrium and to determine the resulting net cooling or heating rates. 

Figure |3] shows examples of net heating/cooling func tions in the presenc e of a photoionizing background. These 
functions were calculated using the MappingsIII code of lAUen et al.l (2008 ) with a photoionizing background com- 



puted self-consistently from the galaxy formation model of iBenson et al.l (l2002ah and an assumed metallicity of 
Z = O.SZq. Unlike coUisional ionization equilibrium cooling functions, which are density independent, photoioniza- 
tion equilibrium cooling curves depend on the density of the gas. The curves shown are therefore computed for gas at 
densities typical of gas in dark matter halos at each redshift indicated. It can be seen that heating becomes important 
for temperatures T <3 x lO^K. 



3.3.2. Heating from Feedback 

The overcooling problem has been a longstanding issue for galaxy formation theory. Simply put, in massive dark 
matter halos simple estimates suggest that gas can cool at a sufficiently large rate that, by the present day, galaxies 
much more massiv e than any observed will have formed ( White and Rees 1978 : White an d Frenk'l99l': 'Katzl 1992 : 
Benson et al.ll2003h . An obvious way to counteract this problem is to heat the cooling gas, thereby offsetting the 



Since the simple scaling is factored out of the cooling function by definition. 
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effects of cooling. Some sort of feedback loop is attractive here in order to couple the heating rate to the cooling rate 
and thereby balance the two. 

In the massive halos where overcooling is a problem energy input from supernovae (see Section refsec:SNeFeedback) 
is insufficient to offset cooling Consequently, much interest has been recently given to the idea that heating caused 
by AGN may be responsible for solving the overcoohng problem . The amount of energy available from AGN is up to 
a factor 20-50 higher than from supernovae dBenson et al ] l2003h . Semi-analytic models have demonstrated that such 



a scenario can work (see ^6.3.2b — producing a significant reduction in the mass of gas able to cool in massive halos 
while simultaneously producing black holes with properties consistent with those observed — under the assumption 
that energy output from AGN can efficiently couple to the coo ling atmosphere in the surrounding dark matter halo 
dCroton et al.ll2006HBower et aLlbOOdlSomerville et al.ll2008bl) . 



The mechanism by which energy output from the AGN is coupled to the surrounding atmosphere remains unclear. 
Several solutions have been proposed, fr om effervescent heating (in which AGN jets inflate bubbles which heat the in- 



tra cluster medium (ICM) via Pd V work: lRovchowdhurv et al.ll2004l:rVecchia et aDl2004l:lRuszkowski et alJl2008l but 
see Vernaleo and Revnoldsl200"6l) to heating by outflows driven by super- Eddington accretion ('Kin22009') and viscous 
dissip ation of sound waves jRuszkowski et alJ2004) with turbulence playing an important role (.Briiggen and Scannapiecol 
200%. 



3.3.3. Preheating 

Heating that occurs prior to the collapse into a dark matter halo can also significantly affect the later evolution 
of baryons. Prior to virial collapse and any associated shock heating the gas evolves approximately adiabatically, 
maintaining a constant entropy. In cosmological studies the entropy is usually written as 

K^ksT/i^mHpl'\ (34) 

where y is a adiabatic index of the gas. An early period of heating (perhaps from the first generation of galaxies) 
can increase the entropy of the gas by raising its temperature. Since entropy is conserved this preheating will be 
"remembered" by the gas. We can consider what happens to such gas when it accretes into a halo, assuming for now 
that a shock forms close to the virial radius. A shock randomizes the ordered infall motion of the gas and therefore 
increases its entropy. At early times, when the shock is relatively weak, this entropy gain will be small compared to the 
preheated entropy of the gas. At late times, the shock entropy will dominate. The result is that the entropy distribution 
of preheated gas in a halo looks similar to that of non-preheated gas, except that there is a floor of minimum entropy. 

The equation of hydrostatic equilibrium can be re- written in terms of entropy using the fact that p = (P/ZT)'^'' 
giving 

dP _ GM(<r)(P\^'y 
dr r 

and shows that the gas will arrange itself in the halo with the lowest entropy material in the center. The presence 
of an entropy floor leads to a density core in the halo center This increases cooling times in the halo core and may 
therefore help prevent the formation o f supermassive galaxies via the overcooling problem (see ^3.3.21 Borgani et alJ 
200lHVoitetal]l2003 : Younger and Brvaa2007.) . Additionally, X-ray observations of galaxy clusters seem to suggest 



\k) ' '''^ 



the presence of entropy floors. 

While preheating is an attractive way of explaining entropy floors in clusters and reducing the overcooling problem 
it is not clear whether it can consistently explain the abundance of galaxies in lower ma ss halos. If preheating oc - 



curred uniformly, it would drastically reduce the number of lower mass galaxies forming (IBenson and Madaull2003l 
Preheating would need to occur preferentially in the sites of proto-clusters to avoid this problem — not inconceivable 
but not demonstrated either. 



'The total binding energy of gas in a dark matter lialo scales as M^'^ while the energy available from supernovae, assuming 100% efficient 
conversion of gas into stars, scales as M. Consequently, in sufficiently massive halos, there will always be insufficient energy from supernovae to 
heat the entire halo. 
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3.3.4. Thermal Conduction 

Shoc k heated gas in dark matter halos is an ionized plasma and is therefore expected to have a thermal conductivity 



given by (lSEitzeflll962D: 



where n^ is the electron density and 



^s = 32(^)"'^rV\ (36) 



Inn^Ze^ In A 

vo = , (37) 



where Z is the atomic number of the ions. If cooling in the inner regions of the hot gas halo leads to a temperature 
gradient in the atmosphere then conduction will cause energy to be deposited at radius r at a rate per unit volume E 
given by: 

d / 2^T\ 



^-^-d^Kd^j' ^''^ 

where T is the temperature in the atmosphere. Thus, thermal conduction can, in principle, act as a heat source for 
the inner regions of the hot atmosphere — as they cool, heat is conducted inwards from the outer, hotter regions of 
the atmosphere. The actual conductivity may be substantially reduced below the Spitzer value if the hot atmosphere 
is threaded by tangled magnetic fields (such that electrons cannot directly transport heat but must effectively diffuse 
along field lines). 

While thermal conduction is an attractive mechanism by which to solve the overcooling problem (since the con- 
ductivity is higher in more massive clusters due to the strong temperature dependence of k^) it ha s been demonstrated 
that it canno t suffici ently offset coo ling rates even with conductivities close to the Spizter value ( Benson et alj 2003 



Dolag et al] 2004; Po pe et al. 2005h. nor can it maintain a stable hot gas atmosphere over cosmological timescales 



jConrov and Ostrikeil 120081: IParrish et aill2009h 



3.4. Fate of Cooling Gas 

In the simplest picture, gas which cools sufficiently below the virial temperature loses pressure support and flows 
smoothly towards the minimum of the gravitational potential well, settling there to form a galaxy. Such a picture is 
lik ely oversimplified, however . 



Mailer and BuUockl (120041) consider the consequences of the thermal instability in cooling atmospheres. They 
find that cooling gas fragments into two phases: cold {T ^ 10^), dense clouds in pressure equilibrium with a hot 
(approximately virial temperature), diffuse component which can persist for cosmological periods of time due to a 
long cooling time. The masses of the clouds are determined from the thermal conduction limit, known as the Field 
lengtl{3 (Field 1965), and processes such as Kelvin-Helmholtz instabilities and conductive evaporation which act 
to destroy clouds. This significantly alters the manner in which fueling of galaxies occurs. The rate of gas supply 
to a forming galaxy now depends on the rate at which the dense gas clouds can infall, due to processes such a 
hydrodynamical drag and cloud-cloud collisions. The timescales for these two processes are given by 

T,™ = — ^ 2.6Gyr (39) 



and 



- 2.4Gyr (40) 



respectively, where md, r^x and Vd are the characteristic mass, radius and velocity of cold clouds respectively, is 
the total mass in clouds, Cd is a drag coefficient, ph is the mean density of hot gas in the halo and is the radius within 
which gas is sufficiently d ense that it has been able to radiate away all of its thermal energy (i.e. the "cooling radius" 



Mailer a nd Bullock"2004'). The numerical values given in each equation are computed for typical cloud properties 



taken from .Mailer and Bullock (.2004) . 



'Above this length scale thermal conduction can damp temperature perturbations in the intracluster medium and prevent cloud formation. 
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Thes e timescales can significan tly exceed the cooling time for gas in galactic scale dark matter halos. Conse- 
quently, iMaller and Bullockl ( 120041) find that considering the formation and infall of such clouds can significantly 
reduce the rate of gas supply to a forming galaxy (by factors of two or so), particularly in more massive halos. This 
picture has recently been confirmed in numerical experiments by Kaufmann et al. ( 2009f) . 



4. Galaxy Interactions 

The original concept of galaxies as "island universes" is, of course, the truth, but not the whole truth. While many 
galaxies do seem to be evolving in isolation there are clear indications that some galaxies are interacting strongly with 
other galaxies or their larger scale environment. We describe these interactions in the remainder of this section. 



4.1. Galaxy Orbits 

A consequence of the hierarchical nature of structure formation in a cold dark matter universe is that dark matter 
halos are built through the merging together of earlier generations of less massive halos. While for a long time 
numerical sim ulations indicated that all trac e of earlier generations of halos was erased during the merging process 
( Katz and White, 1993i; Summers et al.lll995) it was unde rstood on analytical grounds that this was likely a numerical 
artifact rather than a physical result ( Moore et al. 1996ah. Beginning in the late 1990's, N-body sim ulations clearly 
demonstrated that this was indeed the case ( Tormen et al.l 1998 : Moore et al.'l999l Klvpin et al ] ll999h . Unlike earlier 
generations of simulations, they found that halos can persist a s subhalos within larger halos into whic h they merge. 
The current highest resolution simulations of individual halos (IKuhlen et al.ll2008t ISpringel et al.ll2008h show almost 
300,000 subhalo^ and even show multiple levels of subclustering (i.e. subhalos within subhalos within subhalos. . . ). 
Each subhalo may, in principle, have acted as a site of galaxy formation and so may contain a galaxy which becomes 
a satellite in the host potential. 

These subhalos are gravitationally bound to their host halo and, as such, will orbit within it. A subhalo's orbit can 
take it into regions where interactions affect the properties of any galaxy that it may contain. We begin, therefore, by 
considering the orbits of subhalos. 

At the point of merging, which we will define as the time at which the center of mass of a subhalo-to-be first 
crosses the virial radius of its future host halo, we expect the orbital parameters (velocities , energy etc.) to be of order 
unity when expressed in units of the characteristic scales of the host halo. For example, BensonI (2005) shows that 
the radial and tangential velocities of merging subhalos are distributed close to unity when expressed in units of the 
virial velocity of the host halo (Fig.|4]i. This distribution of velocities reflects the influence of the host halo (infall in 
its potential well) but also of the surrounding large scale structure which may have torqued the infalling subhalo. 

Such orbits will typically carry subhalos into the inner regions of halos. Figure |5] shows the distribut ion of orbita l 
pericenters assuming an NFW halo with concentration 10 and the orbital parameter distribution of B ensonI (12005 ). 
Most orbits initially reach to 40% of the virial radius, but a significant tail have orbits which carry them into the inner 
10% of the halo. 



4.2. Gravitational Interactions 
4.2.1. Mergers 

Orbiting subhalos are gravitationally bound to their host hal os and, as such, rarely e ncounter other subhalos at 
velocities resulting in a bound interaction (lAngulo et al. 2008; So merville e t al. 2008b; We tzel et al. Il2008h . 

To cause gravitationally bound interactions between subhalos and their galaxies typically requires a dissipative 
process to reduce their orbital energies. Dynamical friction fulfills this role and tends to drag subhalos down towards 
the center of their host halo, where they may merge with any other galaxy which finds itself there. The classic deriva- 
tion of dynamical friction acceleration from IChandrasekhan (1 1 943b has been used extensively to estimate dynamical 



"This is a lower limit due to the limited resolution of the simulations. The earliest generations of cold dark matter halos may have masses as 
low as IO^'^Mq (depending on the pailicle nature of the dark matter) while state of the art simulations resolve only halos with masses greater than 
around IO^Mq. The ability of even lower mass halos to survive is a subject of much debate ^Berezinskv et al.l200d : izhao et alJlOOTtlGoerdt et al.l 
l2007l;lAngus and ZhaQ.2Q07.:,Elahi et al.„2009i) . 
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Figure 4: Left-hand panel: The joint distribution of radial and tangential velocities of infalling subhalos as they cross the virial radius of their 
host halo. Velocities are expressed in units of the virial velocity of the host halo. Solid contours indicate measurements from a compilation of 
N-body simulations while dashed lines indicate the fitting formula o fiBensonI j2005t) . Contours are drawn at values of dr f IdVgdVr (the normalized 
distribution function) of 0.01, 0.1, 0.5, 1.0 and 1.4. Right-hand panel: The coiTesponding distribution of orbital eccentricities. The distribution 
peaks close to e = 1 (parabolic orbits) and it is apparent that some subhalos are on unbound (e > 1) orbits. Points indicate the distribution measured 
from a compilation of N-body simulations, while th e solid line indicates the distribution found from the same fitting formula as used in the left-hand 
panel. Reproduced, with permission, from lBensonI J2005I) . 




Figure 5: The distribution of or bital pericenters (in units of the virial radius) in an NFW halo having concentration parameter of 10, assuming the 
orbital parameter distribution of lBensor] )2005h . 
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friction timescales within dark matter halos. For example, iLacev and Cold d 1993b applied this formula to estimate 
merging timescales for subhalos in isothermal dark matter halos, finding: 



/(c)Tdyn / r, 

2B(l)lnAU 



where Tdyn = Rv/Vy is the dynamical time of the halo, the virial mass of the orbiting satellite, 

2x J 
B{x) - erf(x) exp(-x ), 



(41) 



(42) 



InA \n{r^V^ IGm^) = ln(My/my) is the Coulomb logarithm (treating the satellite as a point mass), /(e) encapsulates 
the dependence on the orbital parameters through the quantity e - J/JdE) where J is the angular momentum of the 
satellite and JdE) is the angul ar momentum of a circ ular orbit with the same energy, E, as the actual orbit and is 
the radius of that circular orbit. iLacev and Cold d 19931) found that /(e) - e"-^'^ was a good fit to numerical integrations 
of orbits experiencing dynamical friction. 

Besides the fact the dark matter halos are not isothermal, there are a number of other reasons why this simple 
approach is inaccurate: 

1 . Chandrasekhar's derivation assumes an infinite, uniform medium, not a dark matter halo with a radially varying 
density profile; 

2. There is the usual uncertainty in how to define the Coulomb logarithm (e.g. what is the maximum impact 
parameter for interactions between the satellite and dark matter particles); 

3. Subhalos will experience mass loss as they orbit, changing the timescale; 

4. The host halo itself is constantly evolving; 

5. The host halo is non-spherical and has a non-isotropic velocity dispersion. 

While some of these limitations can be overcome (e.g. mass loss can be modeled; (IBenson et al.Ll2002atlTavlor and BabulL 



20041) ; the Coulomb logarithm can be treated as a parameter to be fit to numerical res ults; results exists for dynami 



cal friction in anisotropic velocity distributions; ( Binnevi Il977bt iBenson et al.L 12004 *)) others are more problematic 



Recently, attempts have been made to find empirical formulae which describe the merging timescale. These usually 
begin with an expression similar to the one in eqn. ( 1411 1 but add empirical dependencies on subhalo mass and orbital 



parameters which are constrained to match results from N-body simulations. Results from such studies (IJiang et al 



20081: iBovlan-Kolchin et al.ll2008h show that the simple formula in eqn. (|4TI) tends to underestimate the timescale 
for low mass satellites (probably because it ignores mass loss from such systems) and overestimates the timescale 
for massive satellites (probably due to a failure of several of the assumptions made in this limit). Alternative fitting 
formulae have been derived from these studies. For example. iBovlan-Kolchin et al.l (12008.) find 



^ (My/my)^ 

T^df = TdynA ,^^^^ ^ ^ exp [ce] 



r.(E) 



ln(l + My /my) ' ' ' [ ry 
with A = 0.216,^7 = 1.3, c= 1.9and<i= 1 .0 while [jiTng et al.l ( jooj ) finds 

0^94^0.60 + 0.60 My 1 



T'df - T"dyn- 



2C 



(43) 



(44) 



my ln(l + My /my)' 

with C — 0.43. A comparison of the iBovlan-Kolchin et al ] (l2008h fit, eqn. (141b and measurements from numerical 
simulations is shown in Fig.|6] 

The consequences of merging for the galaxies involved are discussed in 35.2.11 

4.2.2. Tidal Destruction 

An orbiting subhalo and its galaxy will experience tidal forces which may strip away the outer regions or, in 
extreme cases, entirely disrupt the galaxy res ulting in a stellar strea m (as seems to be happening with the Sagittarius 
dwarf galaxy in orbit around the Milky Way; iBelokurov et al.ll2006h . 
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Figure 6: A compa iison of dynamical friction timescales measured from N-body simulations, Tnierge(sim), with the fitting formula of 
iBovlan-Kolchin et al.l (2008), ^merge(fil): is shown by the colored points, with colors coding for orbital circularity: e — 0.33 (red), 0.46 (green), 0.65 
(black), 0.78 (magenta) and 1.0 (yellow). Lines compare rniergc(fit) with the expectation from eqn. 1411 . labeled rniergc(SAM). The dilferent color 
curves correspond to different choi ces of Coulomb logarithms i n eqn. )4U : ln(l + My /my) (black curves) and ^ ln(l + /niy) (blue curves). 
Reproduced, with permission, from lBovlan-Kolchin et aljj2008h . 



D ^ - ^( "'"^"''M + ajV, (45) 



In a rotating frame in which an orbiting satellite instantaneously has zero tangential velocity, the effective tidal 
field felt by the satellite is 

d_/GMh(r)^ 
dr \ 

where Mh(r) is the mass enclosed within radius r in the host halo and co is the instantaneous angular velocity of 
the satellite. An estimate of the radius, r^, in the satellite subhalo/galaxy system beyond which tidal forces become 
important can be made by equating the tidal force to the self-gravity of the subhalo 

= (46) 

Beyond this tidal radius material becomes unbound from the satellite, forming a stream of dark matter and, potentially, 
stars which continue to orbit in the host potential. 

This simple estimate ignores the fact that particles currently residing in the inner regions of a subhalo may have 
orbits which carry them out to larger radii where they may be more easily stripped. As such, the degree of tidal mass 
loss should depend not only on the density profile of the satellite but also on the velocity distribution of the constituent 
particles. Attempts to account for this find that particles in an orbiting satellite that are on prograde orbits are more 
easily s tripped than those on radial orbits which are in turn more easily stripped than those on retrograde orbits 



dRead et a l. 2006). Additionally, some material will be stripped from within the classical tidal radius, as particles 



which contribute to the density inside that radius may be on orbits which carry them beyond it. This can lead to more 
extensive and continuous mass loss as the reduction in the inner potential of the satellite due to this mass loss makes 
it more susceptible to further tidal stripping. (Kampakoalou and Benson.2007i) . 

4.2.3. Harassment 

A less extreme form of tidal interaction arises when tidal forces are not strong enough to actually strip mate- 
rial from a galaxy. The tidal forces can, nevertheless, transfer energy from the orbit to internal motions of stars in 
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the galaxy, effectively heating the galaxy. The gen eric results of such heating are to cause the galaxy to expand 



and to destroy cold , ordered structures such as disks (Moore et al. I ll996bl 1 1 998t Mayer et al.ll200 1 M lGnedinll2003l 



Mastropietro et al.ll2005a.b) . The harassment process works yia tidal shocking in which the stars in a galaxy expe- 
rience a rapidly changing tidal field along its orbit and gain energy in the form of random motions, leading to the 
system expa nding and becoming dynamically hotter. During such tidal shocks, the energy per unit mass of the galaxy 
changes by (lGnedini2003i) 

(AE) = ^lM{r\ (47) 



where (r ) is the mean squared radius of the galaxy and 



/fid - 



EE 



dt 



1 + 



T ' 

^dyn / 



-3/2 



(48) 



where the sums extend over all n peaks in the density field (i.e. the host halo and any other subhalos that it may 
contain) and over all components of the tidal tensor 



d^d) 

dradr/3 ' 



(49) 



where €) is the gravitational potential. Here, t„ is the effective duration of the encounter with peak n and fdyn is the 
dynamical time at the half-mass radi us of the galaxy. The (1 + T^/t^ term describes the trans ition from the 



impulsive to adiabatic shock regimes (iGnedin and Ostrikeilll999t see also lMurah and Weinberglll997allbll3) . 



4.3. Hydrodynamical Interactions 

While the collisionless dark matter is affected only by gravity the baryonic content of galaxies (and their surround- 
ing atmospheres of gas) can be strongly affected by hydrodyamical forces. 



4.3.1. Ram pressure 

The orbital motion of a subhalo through the hot atmosphere of a host halo leads to a large ram pressure. The 
characteristic magnitude of that pressure 



^ ram 



(50) 



can greatly exceed the binding energy per unit volume of both hot gas in subhalos and interstellar medium (ISM) 
gas in their galaxies. As such, ram pressure forces may be expected to quite efficiently remove the hot atmospheres 
of satellite galaxies, a process with several grim aliases including strangulation and starvation, and the ISM of the 
galaxy. 



The first quantification of this process was made bv lGunn and Gotti (119721) who showed that the ram pressure force 
could remove material from a galactic disk if it exceeded the gravitational restoring force per unit area which itself 
cannot exceed 

= 27rG2;*2gas. (51) 

For a disk of mass Md with gas mass Mg having an exponential surfac e density profile fo r both gas and stars with 
scale length r^, the gravitational restoring force per unit area is given by dAbadi et al ] fl999l) : 



T - -, — xe. 

\Tir\ 

d 



(52) 



where x — rfr^ and /q, I\, Kq and K\ are Bessel functions. 

The mass loss caused by this ram pressure can, in many cases, be further enhanced by related effects, such as 
turbulent viscous stripp ing (Nulsen 1982). This initial estimate ha s been revised and calibrated more accurately using 
numerical simulations ( Abadi et^l] 19991; McCarthy et al. 2008). The process of ram pressure stripping has been 
incorporated into some semi-analytic models of galaxy formation dLanzoni et al.l uOOSi ; iFont et al.l 1200 8^ where it 
plays an important role in mediating the transiti on of cluster gala xies from the blue cloud of star-forming galaxies to 
the red sequence of passively evolving galaxies dFont et al.ll2008 ). 
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5. Galactic Structure 



The preceding sections describe the path by which gas achieves sufficiently high density to allow a galaxy to form. 
The next significant question is what structure that galaxy will have. In the following section we will review the 
process of disk and spheroid formation and make quantitative estimates of the sizes of these structures. 



5.1. Disk Formation 

Disks are a common astrophysical phenomenon and galactic disks owe their origins to the same fundamental pro- 
cess as other astrophysical disks: conservation of angular momentum in a system collapsing under gravity eventually 
leading to arrest of the collapse by rotational support. Understanding the physical properties of galactic disks therefore 
requires knowledge of their angular momentum content and the gravitational potential in which they form. 

The angular momentum of the gas which will eventually form a galactic disk arises in the sa me way as that of 
dark matter halos, namely from tidal torques from surrounding large scale structure IHo vld (1 1 949H ' n The magnitude 



of the angular momentum cont ent of the baryonic compo nent of a halo is close to being a scaled version of that of 
the dark matter — for example, van den Bosch et al.l (120021) find that the distributions of spin parameters and angular 
momentum distributions within individual halos for dark matter and non-radiative gas are very similar. They find 
that there is, however, a significant (typically 30°) misalignment between the dark matter and gas angular momenta 
vector^ Numerical hydrodynamical studies show that the angular momenta vectors of galactic d i sks ar e well aligned 



with the angular momenta of the inner regions of dark matter halos ( Bett et al. , 2009 ; Hahn et aL 2010t). 



The distribution of the angular momentum is less well studied but recent investigations (I Sharma and Steinmetz 



20051) using non-radiative hydrodynamical simulations have shown that the differential distribution of specific angular 
momenta, j, is given by 

1 



1 dM 
M"d7 



ma) 



(53) 



where T is the gamma function, M is the total mass of gas, - jtot/a and y'tot is the mean specific angular momentum 
of the gas. ISharma and Steinmeta (l2005b find that the simulated halos have a median value of a - 0.89. The fraction 
of mass with specific angular momentum less than /' is then given by 



f(<j)^r\a^ 



(54) 



where y is the incomplete gamma function. 

We have assumed so far that the angular momentum of the pre-galactic gas is conserved during collapse. This 
may not be precisely true and indeed did not seem to be in earlier generations of hydrodynamical simulations which 
typically found that disk galaxies lost significant fractions of their angular momentum and, as a result, were too 
small (INa varro and Benz 1991; Navarro and White 1994; Navarro et al. 1995). However, more recent simulations 
( Thacker and Couchmanll200lt Isteinmetz and Navarroll2002h do show approximate conservation of angular momen- 
tum (most likely due to the inclusion of effective feedback in these later generations of simulations) and, furthermore, 
conservation of angular momentum leads to disks with sizes comparable to those observed — certainly the gas cannot 
lose a significant fraction of its angular momentum if it is to f orm a disk of comparable size to observed galactic 
disks. Confirmation of these ideas has recently been made by IZavala et al.l (l2008h who convincingly show in N- 
body+hydrodynamical simulations that the particular feedback prescription used can lead to the formation of a disk 
dominated or spheroid dominated galaxy in the same dark matter halo, with stronger feedback leading to a disk 
galaxy. The angular momentum of the mass in the disk galaxy tracks that of the dark matter halo as a whole, grow- 
ing as expected during the linear regime and remaining nearly constant after halo collapse. With weaker feedback 
a spheroid forms instead. Its angular momentum also grows during the linear regime but then declines rapidly (as 
does that of the inner regions of the dark matter halo) due to angular momentum transfer from dense, progenitor blobs 



An alternative view of tliis process, altliougli fundame ntally based upon the s a me pliysics, involves con sidering the angular mome ntum supplied 
to a halo by accreting satellites and has been developed bv lVitvitska et alj i2002h . lMaller et al] i2002h and lBenson and Bowed I20 id) . 

'^These authors note that a significant fraction (5-50%) of the gas in a halo can have a negative angular momentum relative to the total angular 
momentum of the halo, and suggest that this could lead to the formation of low angular momentum spheroids. 
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to the surrounding, dif fuse dark matter Feedback processes are clearly of crucial importance here — for example. 
Wise and Abel ( 2008bl) show that feedback during the formation of protogalaxies at early times can effectively boost 
the angular momentum content of such systems by factors of three to five. This is achieved via forces associated with 
HIT region and SNe blast waves which produce torques much stronger than the usual cosmological torques. 



1983) 



5.7.7. Sizes 

The physics which determines the sizes of galaxy disk s was originally described by Fall and Efstathioul ( 

and later expanded upon by many authors (see, in particular. 'Mo et aL 1998 : Mao et alj|l998l) . 

Once the angular momentum distribution of that part of the baryonic component of the halo which cools to form the 

galaxy is known, finding the structure of the resulting rotationally supported disk is reduced to solving the following 

equation 

4^ = UiRl (55) 
R\M) dR 

where y(M) is the specific angular momentum enclosing mass M and this equation is solved for R{M). The potential 
is a sum of the self-gravity of the disk and that of any external potential (dark matter halo and bulge for example) 
wh ich may have responded adia batically to the formation of the disk. For a thin disk, the potential is generically given 
by dBinnev and Tremainell2008h 



<D(7?,0) = -4G f , f 
Jo ^/W^ J a 



dR' 



R'I.(R') 



(56) 



which, for an exponential surface density profilJ^ 



Soexp|- — 



simplifies to 



O(/?,0) = 7TG^oR[Ioiy)Ki(y) - Ii{y)Ko(y)], 



(57) 



(58) 



where y - R/IR^ and /o, /i, K^) and Ki are Bessel functions. Knowing the functional form of /'(M) it is possible to 
solve for R(M) and therefore the density profile of the disk. Often, for simplicity, a particular form for the density 
(e.g. exponential) is assumed which leaves a single free parameter (the scale length) to be solved for. 

A significant complication to this picture arises from the fact that the external potential in which the disk forms 
is likely to change in response to the formation of the disk. If the disk forms slowly, such that the timescale for 
changes in the pote ntial greatly exceeds the dy namical time of the dark matter halo, then we can use adiabatically 
invariant quantities dBinnev and Tremainel2008l) to esti mate the response of the h alo to the forming disk. The original 
formulation of this argument for galaxies was given bv lBlumenthal et al.l (119861) and has been used extensively ever 
since. In this simplified picture, the dark matter particles are considered to be on circular orbits in a spherical potential. 
In this case, the only adiabatic invariant whose corresponding angle variable has a non-infinite period is the magnitude 
of the angular momentum, L^. Prior to the formation of any galaxy, this angular momentum is given by 



(59) 



where Mdm('") is the mass of dark matter within radius r, ro is the initial radius of the circular orbit and we have 
assumed that baryons are present in the halo at the universal fraction, fb(= Qb/[Qb + ^^cdm]) and distributed as the 
dark matter If is conserved during the formation of a galaxy, then, at any point after galaxy formation has begun 
we have 



4 = (1 - fO GlM'^^irf) + Mg,i(rf)]rf, 



(60) 



''^These equations assume a razor thin disk. For a thickened disk a density distribution p(R, z) = po exp(-S/Sd)sech^(2/22o), with zo being the 
scale height of the disk, is often assumed as it is the self-consistent d istribution for an isothermal population of stars I Spitzer 1942) . The resulting 
gravitational potential can be found (e.g. iKuiiken and Gilmorel 19891) although the calculation is significantly more involved. 
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where M'(r) is the mass of dark matter enclosed within radius r after responding to the forming galaxy, Mgai(r) is 
the mass of the galaxy within that same radiu^^ is the baryon fraction remaining uncondensed into the galactic 
phase and rf is the radius of the new circular orbit for this dark matter particle. If we assume no shell crossing, then 
Mp^(rf) - Mdm('"o)- Given a knowledge of the mass and density distribution of the galaxy this allows us to solve for 
rf as a function of the initial radius ro. 

This approach makes several simplifying assumptions, including adiabatic growth, a spherical halo and circular 
orbits for dark matter particles. Many of these assumptions have yet to be adequately tested. However, recent work 
has attempted to address the last of these. With a more realistic distribution function for the dark matter, the particles 
will occupy a range of orbital eccentricities. In a spherical potential, the angular momentum is joined by a second 
adiabatic invariant, the radial action: 



L, 



^ Jr 

^ ' mm 



(61) 



where the integral is taken along the orbit of a particle from its pericenter to apocenter The eff ects of taking i nto ac - 
count this second in variant and a physically reasonable orbital distribution have been explored by Gnedin et al. ( 20041) . 
Gnedin et al.l (120041) find that the simple model described above systematically overpredicts the degree of contraction 
of the halo. However, they find that a simple modification in which the combination M(r)r in the above is replaced 
with M(r)r, where ? is the orbit-averag ed radius, approxirn ately accounts for the orbital eccentricities and gives good 
agreement with numerical simulations. iGnedin et al.l (120041) further find that the orbit-averaged radius in cosmological 
halos ca n be well fit by the rela tion r/ryir - ACr/rvii)" with A - 0.85 + 0.05 and w = 0.8 ± 0.02. More recent sim- 
ulations (iGustafsson et al. 20061) support this picture, but show that the optimal values of A and w vary from galaxy 
to galaxy and seem to depend on the formation history (i.e. on the adopted rules for star formation and feedback; see 
also Abadi et al. 2009). This suggests that the process of halo response to galaxy growth is more complicated than is 
captured by these simple models, perhaps because the assumed invariants are not precisely invariant in non-spherical 
potentials or perhaps because galaxy growth is not sufliciently slow to be truly adiabatic. 

Galaxy disks are, of course, not razor thin. The origins of the vertical extent of galaxy disks remain a topic of 
active research but various possibilities are currently considered seriously: 

1. External origins (accretion): 

(a) The hierarchical nature of galaxy formation implies that a galactic disk can expect to accrete pre-existing, 
smaller stellar systems frequently during its life. The stars from these smaller galaxies are often found 
in numeric al simulations to form a thickened disk structur e, in the same plane as the pre-existing disk of 
the galaxy dAbadi et al.ll2003t IVillalobos and Helmill2008b . This is to be expected as dynamical friction 
will tend to drag orbiting satellites into the plane of the disk where tidal forces will proceed to shred the 
satellite, leaving its stars orbiting in the plane of the disk. 

(b) Early, chaotic gas accretion: In a slightly diff'erent scenario. iBrook et al. I (l2004 find, in numerical simula- 
tions of a forming galaxy, that many of the thick disk stars form from accreted gaseous systems during an 
early, chaotic period of merging. 

2. Internal origins (dynamical heating): 

(a) Dark matter substructure: Buildin g upon earlier works (iToth and OstrikeJl992l: l lOuinn et al .11993^ ISellwood et al 
19981: 1 Velazquez and Whitelll999h . in recent years there have been numerous studies focused on the ques- 



tion of whether relatively thin galactic disks can survive in the rather violent environment of a hierarchi 
cally formed dark matter halo ( Font et al. I2OOII; Benson et al. 2004 ; Kazantzidis et al. 2008 ; Read et al 



2008; Kazantzidis et al 1 120091 1 iPurcell et al J l2009t) . A careful treatment of dark matter substructure orbital 



evolution and the build-up of the substructure population over time is required to address this problem. 
The consensus conclusion from these studies is that galactic disks can survive in the currently accepted 
cold dark matter cosmogony, but that interactions with orbiting dark matter substructure must be a signif- 
icant contribution to the thickening of galaxy disks. The heating is dominated by the most massive dark 
matter substructures and so is a rather stochastic process, depending on the details of the merger history 
of each galaxy's dark matter halo and the orbital properties of those massive substructures. 



'Note that we are approximating the galaxy as a spherical mass distribution here. 
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(b) Molecular clouds: Massive molecular clouds can gravitationally scatter stars that happen to pass by them, 
effectively transforming some of their orbital energy in th e plane of the galaxy disk into motion perpendi c- 
ular to that plane, thereby effectively thickening the disk ( Spitzer and Schwarzschildlll953 : Lacevlll984 ). 

(c) Spiral arms: Spiral density waves can also act as scatterers of stars, but predo minantly increase the stellar 
velocity dispersion in the plane of the galaxy, resulting in little thickening dJenkins and Binnev 1 19901: 



Minchev and Ouillenll2006h . 



Most l ikely, a combination of these processes is at work. A recent review of these mechanisms is given bv lWen and Zhao 
(|2004 . 



5.7.2. Stability 

Spiral arms and other non-axisymmetric features such as bars in galactic disks are a visually impressive reminder 
that these system s possess interestin g dynamic s. The basic theory of spiral a r ms, that they a re density waves, was 
first proposed by iLin and Shul (Il964i see also iMarochnik and Suchko vl l 19961: 



Binnev and Tremaine 2008i) and has 



come to be widely accepted. From our current standpoint, the question in which we are interested in is whether these 
perturbations to an otherwise smooth disk are stable or unstable and, if unstable, how they affect the evolution of the 
gal axy as a whole. 



Toomrd (Il964t see also iGoldreich and Lvnden- Beiil l 1965h derived an expression for the local stability of thin 



disks to axisymmetric modes in the tight-winding approximation which turns out to be extremely useful (and often 
approximately correct even in regimes where its assumptions fail). Disks will be unstable to axisymmetric modes if 
Q < 1 where 



Q 



where 



dQ^ 



gas 



•4Q' 



1/2 



(62) 



(63) 



is the epicyclic frequency and D. is the angular frequency of the disk, for a gaseous disk of surface density Zgas and 
velocity dispersion cTgas, and 



Q 



(64) 



3.36G5:* 

for a stellar disk of surface density and velocity dispe rsion cr^r. For a disk co n sisting of two com ponents, 
gas and stars, a joint stability analysis was carried out by Jog and SolomonI ( 1984 ). Efstathiou ^200&. see also 
Bertin and Romeo, 1988.) solves the resulting cubic equation for the most unstable mode and finds a criterion 



' gas 



(65) 



where a and /3 are t he ratio of stellar t o gas velocity dispersio ns and surface densiti es respectively and the function 
g{a,/3) is as given by Efstathiou ( 2000h and is shown in Fig.lTl IWang and SiUd ( 19941) also show that the stability of a 
two-component disk can be approximated by 



1 \-i 



(66) 



while a more general criterion taking into account the thickness of the disk was found bv lRomed (119921) and lRomeo 
(11994'). 

Toomre's criterion applies to local perturbations. Perturbations on the scale of the disk can occur also. Study of 
these global instabilities is less amenable to analytic treatment (since one can no longer ignore contributions from 
distant parts of the disk). Significa nt work on this subject began in the 1970's. The classic result from that time is 
due to Ostri ker and PeeblesI (Il973h who found that rapidly rotating, self-gravitating stellar systems would become 
violently unstable to non-axisymmetric m - 2 modes if T/\W\ >0.14 (for a stellar system) or 0.27 (for a fluid system) 
where T and W are the kinetic and potential energies of the system. Such systems can be unstable to global pertur- 
bations even if they satisfy Toomre's stability crite rion. These values were app roximately confirmed by numerous 
numerical studies, but also challenged by others (see lChristodoulou et al. (1995b) for a discussion of this). Later work 
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Figure 7: The function g{a,fi) appearing in eqn. (65). g{a,fi) is shown as a function of /} (the ratio of stellar to gas velocity dispersions) for various 
values of a (the ratio of stellar to gas surface density). Reproduced, with permission, from Efs tathiou ( 2000). 



( Christodoulou et al. I ll995ah lias demonstrated that r/|lV| is not a sufficiently general parameter to determine stabil- 
ity in a ll systems. More recent studies have proposed improved stability criteria. For example IChristodoulou et al 
(Il995bl) find that 



T/\W\ 



(67) 



where Qj is the Jeans frequency in the radial direction and Q. is the mean angular velocity of the system, is a much 
better indicator of stability. This can be approximated by 



a = Tj/\W\, 



(68) 



where Tj - LQ.j/2 and L is the total angular momentum of the system. IChristodoulou et al. (Il995bh found that 



a < 0.254-0.258 is required for stability in stellar systems while a < 0.341-0.354 is required for gaseous systems. 
IChr istodoulou et al. (1995b) also demonstrate that this criterion is approximately equivalent to the alternative form 
found bv .Efstathiou et al.. (.1982 ) on the basis of numerical simulations of exponential stellar disks: 



>1.1, 



(69) 



and show that an equivale nt result for gaseo us disks would replace the 1.1 with 0.9. A much more extensive review 
of disk stability is given bv lSellwoodl(l2010l) . 

Global instabilities most likely lead to the formation of a very strong bar which effectively disrupts the disk 
leaving, after a few dynamical times, a boxy /peanut b ulge or a disky bulge (A thanassoula 2008). This may therefore 
be a possible formation scenario for pseudo-bulges dKormendv and Kennicutt) 12004 ). i.e. bulges formed through 
secular processes in the disk (see ^5.2.2b rather than as the result of a merger event (see ^5.2. lb . 



5.1.3. Bars/Spiral Arms 

The presence of perturbations in the disk can also affect its structure in less dramatic but still significant ways. 
(Once again, the work by Sellwood (2010) gives a more in dept h review of these subjec ts and is highly recommended.) 
For example, numerical simulations by iRoskar et aP (l2008bt see also iRoskar et al ] l2008al) show that a significant 
fraction of stars in galaxy disks undergo large migrations in radius (as shown in Fig.|9]l due to resonant scattering from 
spiral arms. The resonant scattering of a star initially on an approximately circular orbit tends to change its energy 
and angular momentum in such a way as to keep it in a circular orbit, but to change the radius of that orbit. In a frame 
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Figure 8: The Lindblad diagram showing allowed regions of the energy-angular momentum plane (white, with forbidden regions shaded gray) 
bounded by the curve corresponding to circular orbits. Arrows indicate how scattering from spiral density waves moves stars in this diagram. At 
the corotation resonance (CR) this motion parallels the circular orbit curve, implying that scat tering at corotat i on of stars initially on circular' orbits 
leaves those stars on approximately circular orbits. Reproduced, with permission, from .Sellwood and BinnevI i2002h . 

cor otating with the spiral pertur bation there is an energy- like invariant quantity know as Jacobi's integral, Ej, defined 
as dBinnev and Tremainell2008h 

Ej^E- QpL, (70) 

where E and L are the energy and angular momentum of an orbit respectively and Op is the pattern speed of the 
perturbation. This implies that changes in energy and angular momentum are related by 

AE = QpAL. (71) 

As shown in Fig.|8]the slope of the E-L relation for circular orbits at the corotation resonance is f2p. As such, stars 
initially on c i rcular o rbits near the corotation resonance will tend to be scattered into approximately circular orbits. 



Roskar et al. show that this phenomenon can lead to the formation of an outer stellar disk with a steeper 



exponential decline than the inner disk and formed entirely of stars scattered from the inner regions of the galaxy. 
This process also smooths out age and metallicity gradients in disks that otherwise form inside out and so tend to be 
older and less metal rich in the center. 

5.2. Spheroid Formation 

The formation of galactic spheroids (which we take here to mean both the bulges of disk galaxies and isolated 
elliptical galaxies) can proceed via two distinct routes. The first, through the destruction of pre-existing stellar systems 
in violent mergers, is a natural consequence of hierarchical galaxy formation. The second, secular evolution of galactic 
disks, is a natural consequence of the dynamics of self-gravitating disk systems. It is worth noting th at, observationally, 
it has been proposed that bulges be divided into two broad classes dKormendv and Kennicuttl200 4*): "classical" bulges 
(those which look like ellipticals in terms of their light distribution and kinematics but happen to live inside a disk) 
and "pseudo-bulges" (those which do not look like elliptical^^)- There is good evidence (see ^4.2.1l i that ellipticals 
(and therefore presumably classical bulges also) form through major mergers, and suggestions that pseudo-bulges are 
the result of secular evolution of galactic disks. 



'^Typically, pseudo-bulges are better fit by a lSersid j 19681) profile with index « « 1 (as opposed to the n = 4 profiles of classical bulges), have 
significant rotation and may show signs of disk phenomena such as bars or spiral features. 
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Figure 9: The distribution of formation radii, Rfonn (in kpc), for stars in the Solar neighborhood. The black histogram shows the distribution from 
all st ars while the red and blue histograms show the contributions from high and low metallicity stars respectively. Reproduced, with permission, 
from lRoskar et alj <2008d) . 



5.2.1. Major Mergers 

As was discussed in ^4.2. II dissipative processes can lead to gravitationally bound interactions between galaxies. 
When the masses of the galaxies in question are comparable we may expect significant changes in their structure and 
the formation of a merger remnant which is very different from either of the merging galaxies. Such major mergers 
are thought to be responsible for the formation of spheroidal (i.e. elliptical) galaxies from pre-existing galaxies (which 
could be, in principle, disk galaxies, other spheroidals or somethi ng intermediate). 

The process of violent relaxation ( Lvnden-Belll 1 9671: see also Tremaine et al. ill 9861) . in which the energy of orbits 
undergoes order unity changes due to signficant time-variable fluctuations in the gravitational potential, leads to a 
randomization of the orbits leading to a Maxwellian distribution of energies but with temperature proportional to 
stellar mass. This can turn the ordered motions of disks into the random motions seen in spheroids. This process is, 
however, rather poorly understood — it seeks an equilibrium state which maximizes the entropy of the system, but the 
usual entropy is unbounded in gravitating systems implying no equilibrium state exists. Arad and Lvnden-Belll ( 2005 ) 
demonstrate this problem by showing hysteresis effects in violently relaxed systems (i.e. the final state depends on 
how the system goes from the initial to final states). 

The remnants of major mergers of pu rely stellar disk systems, while spheroidal, do not look like elliptical galax- 
ies. As shown bv iHernquist et al.l (Il993h their phase space densities are too low in the central regions compared to 
observed ellipticals. This implies that mergers between reasonably gas rich (gas fractions of around 25%) galax- 
ies are required— the pr esence of gas allows for dissipation and the formation of higher phase-space density cores. 
Robertson et al. 120069) find that similar gas fractions and the subsequent dissipation are required to produce the 



observed tilt in the fundamental plane of elliptical galaxies — mergers of purely stellar systems instead follow the 
expected virial scalings. 

Mergers are often separated into major (mergers between galaxies of comparable mass) and minor (mergers in 



which one galaxy is significantly less massive than the other). Numerical simulations (e.g. Bour naud et al 



2005) 



show that mergers with a mass ratio ju = M2/M1 >0.25 are able to destroy any disks in the ingoing galaxies and leave 
a spheroidal remnant, while mergers with lower mass ratio tend to leave disks in place (although somewhat thickened). 
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As this transition is unlikely to be a sharp one. lSomerville et al. I (l2008bh propose a more gradual transitiorFl with 
fraction 



/sph - 1 - 



1 + 



A* 



ellip 



(72) 



where /eiup ~ 0.25-0.30, of the disk stellar mass being put into the spheroid component of the merger remnant. 

Mergers are expected to trigger a, possibly very large, enhancement in the star formation rate in th e merging 



system . The strengths of such bursts of star formation depend upon the properties of the merging galaxies. ICox et al 



(l2008h have calibrated this enhancement against a large suite of N-body+hydrodynamics merger simulations. They 
define the "burst efficiency", e, to be the fraction of gas consumed during the interaction above that which would have 
been consumed by the constituent galaxies in isolation during the same time period. They find that the burst efficiency 
is well fit bxTl 



\ jr 





if Msat/Mprimary > ^0 
if Msat/Mprin,ary < 



(73) 



where ci-a = 0.56, y = 0.50 and eo = 0.09. 

It has recently become clear that not all major mergers lead to the formation of a spheroid. Under certain co ndi- 
tions, major mergers of very gas rich systems can lead to th e reformation of a disk after the merger is over (Bar nesI 
20021: IStTi-ingel and Hernquistll2005l: iRobertson et alJ l2006ah . This requires a high gas fraction (greater than about 
50%) just prior to the final coalescence of the merging galaxies and therefore may preferentially occur under condi- 
tions which prevent the rapid depletion of gas after the first passage of the galaxies. 

5.2.2. Secular Evolution 

Major mergers are not the only w ay to form a spheroid. Internal, secular processed in galaxies can also disrupt 



the cold and relatively fragile disks (IKormendv and Kennicutt l2004 . In particular, bars (a disk phenomenon) can 
efficiently redistribute mass and angular momentum and lead to the build-up of dense central mass concentrations, 
reminiscent in many ways of bulges formed through mergers. To distinguish such se cularly formed bulges from thei r 
merger-formed (or "classical") counterparts, they are referred to as "pseudo-bulges" ( Kormendv and Kennicuttl2()04l) . 

Such secular processes are the result of quite generic dynamical consideration^ and so most likely operate in all 
galaxies. Whether or not they are important depends upon their timescale. For example, relaxation due to star-star 
encounters in a galaxy operates on a timescale many orders of magnitude longer than the age of the Universe and so 
can be safely neglected. Instead, most relevant secular processes involve the interaction of stars (or gas elements) with 
collective phenomena such as bars or spiral arms. 

A general picture of how secular evolution leads to the formation of pseudo-bulges has emerged. As a bar spon- 
taneously begins to fornP*! it transfers angular momentum to the outer disk and increases in amplitude. The response 
of gas to this bar is crucial — gas accelerates as it enters the bar potential and decelerates as it leaves. This leads to 
shocks forming in the gas which lie approximately along the ridge line of the bar. These shocks lead to dissipation of 
orbital energy and, consequently, inflow of the gas. The enhancement in the gas density as it is concentrated towards 



' ISomerville et all t2008bh g ive no justification for tliis particular functional form, it is merely intended to give a smooth but rapid transition. 

'^The equation in lCox et id is missing the coefficient of en (Cox, private communication). 

"Generically, any internal dynamical process operating on a timescale significantly longer than the dynamical time. 

^"As pointed out bv l Kormendv and Kennicutl i2004) disks are fundamentally prone to spreading in the presence of any dissipative process, where 
mass is tra nsported inwards and angular momentum outwards, because this lowers the energy of the system while conserving angular momentum 
iLvnden-B ell and Kalnais 1972). This result c an be traced back to the negative specific heat of gravi tating systems, and is analogous to the process 
of core collapse in three dimensional systems iLvnden-Bell and Kalnaisll97 2: Binnev and Tremaine 2008). 

^'The bar instability in galactic disks involves some fascinating dynamics. Binnev and Tremaine (2008) give a clear explanation of the physics 
involved. Briefly, the bar instability involves the joint actions of the swing amplifier and a feedback mechanism. A randomly occurring leading 
spiral density wave in a disk will unwind and, as it does so, will rotate faster. As it swings from leading to trailing it reaches a maximum rotation 
speed which is close to the average orbital speed of stars in the disk. This leads to a resonance condition, in which the wave can strongly perturb 
the orbits of those stars, the self-gravity of which enhances the bar further, leading to an ampHfication of the wave. If there is some mechanism to 
convert the amplified trailing wave that results into a leading wave once more (e.g. if the wave can pass through the center of the galaxy and emerge 
as a leading wave, or if nonlinear couplings of waves can generate leading waves) the whole process can repeat and the wave will grow stronger 
and stronger. 
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the galactic center inevitably leads to star formation and the build-up of a pseudo-bulge. Bars eventually destroy 
themselves in this way — the increase in the central mass of the galaxy effectively prevents the bar instability from 
working. 

While this general picture seems to be incontrovertible, many of the details (e.g. the redshift evolution of bars, 
the importanc e of the dark matter halo and the cosmological evolution of the disk) re main poorl y understood, but ac- 
tively studied (Friedli and BenJl993l; Shlosman and Noguchiil993;,El-Zant and Shlosmaa2002 : Athanassoulal2002 : 



Bournaud and Combesll2002l;lshen and Sellwoodll2004 



5.2.3. Sizes 

The sizes of spheroidal galaxies formed through major mergers can be determined given the properties of the 
progenitor galaxies and their orbit and some knowledge of the extent to which mass and energy is conserved through 
the merging process. If the galaxies are dissipationless, purely stellar systems (a so-called dry merger) then energy is 
at least approximately conservecQ In a gas rich (or wet) merger the gaseous component is likely to radiate significant 
amounts of energy prior to forming stars. 

A simple model to compute the size of a galaxy formed via a major merger was described bv lCole et al. 
who assumed perfect conservation of mass and energy. The size of the merger remnant, rf, is then given by 



CfG(Mi + M2)2 
rf 



C2GMI 



+ 



n r2 
/o rbitGMiMa 
ri+ r2 



(74) 



The left-hand side represents the gravitational binding energy of the final system while the first two terms on the right 
represent the gravitational binding energy of the two merging galaxies, which have masses Mi and M2 respectively 
and half-mass radii r-i and r2 respectively. The "c" coefficients relate the actual bind i ng ene rgy to the characteristic 



value of GM^/r and depend on the density distribution of the galax}]^ ICole et aL find c = 0.45 for a de 



Vaucouler's spheroid and c — 0.49 for an exponential disk and so adopt c - 0.5 for all galaxies (which may be 
composite systems consisting of both a disk and a spheroid) for simplicity. The final term in eqn. ( l74l i represents the 
orbital energy of the two merging galaxies just prior to merging and is parameterized in terms of their gravitational 
binding energy when separated by the sum of the half-mass radii. Cole et al. (2000) adopt /orbit = 1 so that this orbital 
term corresponds to a circular orbit. According to this approach, a merger of two identical galaxies of radii ri results 
in a remnant of siz e rf = (4/3)ri. 



More recently, ICovington et al. 1 (l2008h performed numerical simulat ions of me r ging g alaxies including gaseous 
components, and found significant deviations from the simple model of Cole et al. ( 2000h as a result of energy loss 



due to radiative processes in merging galaxies containing gas. Covington et al.l ( 20081) propose an improved model 
which accounts for this energy loss and results in a merger remnant size given by 



^int.f ~ ^int.i "t" ^rad "t" 



(75) 



where 



is the internal binding energy of the final galaxy. 



int.f 



cGM^ 



rf 



cGM^ 
ri 



cGMl 
r2 



(76) 



(77) 



In principle, energy may be lost from the stellar system by being transferred to the dark matter, or high energy stars may be flung out beyond 

the escape velocity, removing both mass and energy from the system. 

^^^In principle, the binding energy terms should also account for the stellar-dark matter binding energy. Cole et all jjOOff) ignore the dark matter 
contribution of the binding energy, effectively assuming that its contribution to each term in eqn. I74t scales in proportion to the stellar masses of 
the galaxies. 
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is the internal binding energy of the initial galaxies. 



GM1M2 1 , 1 , 

^ Sep 



(78) 



is the orbital energy of the galaxies just prior to the merger when their separation is rjep and their center of mass 
velocities are Vi and V2, and 



(79) 



is the radiative energy loss where Ki and fgj are the initial internal kinetic energy and gas fraction for progenitor 
The impulse, ,■ is defined to be = AE/Kto^ where K^ot is the total initial internal kinetic energy of the galaxy and 



AE = 



.(r^ . + tot'-peri + Cr?, ,) ' 

pen^ pen ^'^"^ l,tot^ 



(80) 



where Mi.tot and Ma.tot are the total masses of the galaxies (baryonic plus dark matter), ri tot is the total half-mass 
radius of the primary galaxy and rperi is t he pericentric distance of the first passage of the two galaxies, treating them 
as point masses from their initial orbit. ICovington et alj (120081) find that values A = 1 .6, B = 1 .0, C = 0.006 and 
Crad = 1 best fit their simulation results. 



Similar arguments have been applied to the formation of spheroids through secular processes (e.g. ICole et al 

but the assumptions of this approach are far less well tested in such cases. 



6. Star Formation, AGN and Feedback 

So far, we have not considered how a galaxy converts its available gas into stars and what effect those stars, and 
associated supernovae explosions, may have on the evolution of the galaxy. Additionally, observational evidence 
suggests that all galaxies contain a supermassive black hole at their center, which may play a key role in limiting the 
process of galaxy formation. We will discuss these aspects of galaxy formation theory b elow. We do not expl ore in 



detail the role of the associated galactic winds in galaxy evolution, but refer the reader to lVeUleux et alJ (l2005h for a 
thorough review. 

6.1. Star Formation 

Star formation theory warrants an entire review to itself (for a recent example see McKee and Ostrike3l2007 ) and 



so we will summarize only those aspects most pertinent to galaxy formation theory. The past decade has seen a 
greatly improved understanding of how the key processes of turbulence, magnetic fields and self-gravity interact to 
form molecular clouds and stars. Nevertheless, there remain numerous unsolved problems in star formation theory. 
These problems propagate into galaxy formation theory if we wish to understand the rate at which stars form in a 
galaxy and any consequences that may have for further galactic evolution. This problem is somewhat mitigated by 
the fact that, for galaxy formation theory, we do not necessarily care about the details of how stars form. Instead, we 
would simply like to know, given the large scale properties of a galaxy (e.g. mass, size, density, dynamical time, gas 
fraction, chemical composition), what is the resulting rate of star formation. 

Traditionally, this question has been answered by appealing to empirical r ules or d i mensi onal analysis. For ex- 



ample , much use has been made of the empirically derived scaling relations of' Schmidt (Il959h and lKennicuttl (Il989[ 



19981) . The Schmidt-Kennicutt law states that the rate of star formation per unit surface area, 2*, depends on the 



surface density of gas, Egas as 

U^I-l,,, (81) 

where n is measured to be approximately 1 .4. While this relation is practically useful for galaxy formation theory, 
insofar as it allows one to bypass the question of star formation and move directly to a star formation rate, it suffers 
from the same issue as all empirical relations used in theoretical models: there is no way to be sure that it is valid 
beyond the regimes where it was originally measured. 
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Recently, new observations have suggested a relatively simple model for star formation dKrumholz and Tanl2007h . 
Stars form in molecular clouds and so it is natural that it would be the density of molecular gas (rather than total) gas 
which matters. Furthermore, it is observed that a constant fraction (around 2%) of molecular gas turns into stars per 
free-fall timescale. This implies - eff/HjPgas/fff where eff is the fraction of molecular gas turned into stars per 
free-fall time, /h, i s the molecular fraction and % i s the loca l free-fall timescale. 

In the model of Robertson and Kravtsov ( 2008 , see also iKrumholz et alj|2"009l) the star formation rate is given by 



2* MhsFRh lf'^ll, 



(82) 



where /h, is the fraction of hydrogen in mole cular form, /isfr is the s cale height of the star-forming gas and hg^s is 
the scale height of the interstellar medium. Robertson and Kravtsovl (12008 ) compute /h^ and the scale height self- 
consistently in the presence of the interstellar radiation field produced by earlier generations of stars. This results in a 
net scaling with total gas surface density of 

^ (83) 

where ntot ~ 2 for massive galaxies and ntot >4 for dwarf galaxies. When expressed in terms of the molecular hydrogen 
gas surface density a scaling 

i*ocE^"»', (84) 

with Hmoi ~ 1.3 is found for all galaxies. These scalings are consistent with a broad range of observations of star 
formation in d ifferent galaxies . 

Similarly, IKrumholz et al. I (l2009h find that the following expression encapsulates the complex physics of star 
formation: 



2* 



/H3(2g, c,Z') 



2.6 Gyr 

0.33 



( 85MgPC-2 ) 



where the molecular hydrogen fraction is 

/H,(2g,c,Z')^ 1 



85MopC- 
85MopC- 



< 1 
> 1, 



(85) 



4 1-1-5 



-1/5 



(86) 



and* = ln(l+0.6;(f)/(0.042;eomp,oZ'),A' = 0.77(l+3.1Z'"-^''5)^ ^ ^ 0.07 12(0. li-'+0.675)-^Meomp,o = Xcomp/(lMopc-^) 
and Z' is the metallicity normalized to the Solar value. The transition at gas surface densities of 85Mopc"^ corresponds 
to the point at which the ambient pressure becomes comparable to the internal pressure in molecular clouds (and 
therefore influences the properties of those molecular clouds). Here, Scomp is the surface density of a ~ 100 pc-sized 
atomic-molecular complex. Since simulations and semi-analytic models typically only predict the gas distribution on 
scales significantly larger than this we can write Scomp - cSg where c is a clumping factor that accounts for structures 
which are unresolved in the simulation. This clumping factor should approach unity as the resolution approaches 100 
pc at which point molecular cloud complexes should be adequately resolved. 



6.2. Black Hole Formation 

Over the past ten years it has become possible to measure the masses of supermassive black holes residing at 
the centers of galaxies for relatively large samples. The existence of strong cor relations between the m asses of 
these black holes and the properties of their host galaxy — such as spheroid mass ("Magorrian et al."l998'), velocity 
dispe rsion (Ferrarese and Merritt 20 00; Gebhardt et ^ 1, 2000a,b), number of globular clusters (Burkert an d Tr emaind 
20101) or even host dark matter halo (lFerraresell2002l) — is suggestive of som e interaction be twe en forming galaxies 
and su permassive black holes. Of course, correlation does not imply causation ('Munroe'2009') and'jahnke and Ma cciol 
(l2010h show that a black hole-host galaxy mass relation can arise from uncorrected initial conditions via simple 
merging, but the theoretical need for large amounts of energy to inhibit galaxy formation in massive halos naturally 
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leads to the idea that supermassive black holes (SMBH) and galaxy formation are connected dBenson et alj|2003h . 
In addition, understanding the formation of these most massive of black holes is interesting in its own right and 
has important observational consequences for both studies of active galactic nuclei and gravitational wave detection 
experiments such as the Laser Interferometer Space Antenn^\ 



galaxies ( Silk and Rees»1998;.Monaco et albOOa.Kauffmann and Haehneltl200Q; Kine 2003;.Wvithe and Loebl2003; 


Matteo et al. 2003; Volonteri et al. 2003 


; Kind 2005; 


Beaelman and NathI l2005t Malbon et alj 20071 Matteo et al. 


20081 Volonteri et all 20081 Siiacki et al. 


2009 ; 1 Volonteri and Nataraianll2009h. We will not review all details of this 



process here. Instead, we will focus on the basic formation mechanisms and, in 36.3.21 on the interaction between 
SMBHs and galaxy. 

Before supermassive black holes can grow via accretion or merging, there must be some pre-existing (probably 
not supermassive) seed black holes. Most plausibly, these seeds form at high redshifts as the remnants of the earliest 
generation of Population III stars which have reached the end of their stellar lifetimes. Details of the formation of 
these first stars remain incompletely und erstood, but hydrodynamical si mulations suggest that they have masses in 
the range of a few hundred Solar masses dAbel et al.ll2002 : Bromm et al.lE002) . leaving intermediate mass black hole 
remnants. 

To determine the rate at which gas accretes onto a black hole (or, more precisely, the black hole and any associated 
accretion disk system) we must consider how the black hole afifects the gas through which it is moving. The gravity 
of a fast moving black hole will deflect gas that passes by it, focusing it into a wake behind the black hole which will 
then accrete on to the black hole. Th is problem, and the equivalen t for a slowly (subsonically) moving black hole, was 
first studied bv lHovle and Lvttletonl(il939i) and lBondi and Hovld (1 1944)) . This leads to an accretion rate of 



Mr 



(c2 + v2)3/2' 



(87) 



where M, is the mass of the black hole, Cj is the sound speed and v the relative velocity of black hole and gas. The 
accretion occurs from a characteristic radius of 



'"BHL 



GM. 



(88) 



An in-depth review of Bondi-Hoyle-Lyttleton accretion is given bv lEdgad (120041) . The growth of black holes will 
be enhanced by any process which increases the density in the central regions of the galaxy in which they reside. 
At early times, this ma y occur due to gravitationally unstable disks forming bars and driving gas towards the cen- 



ter (IVolonteri et al.ll200 8). while at later times galaxy-galaxy mergers can results in dissipation and gas flows to the 



center jMatteo et al.'' 2005|). O ther mechanisms for delivering gas to the centers of galaxies are also possible, for ex- 
ample|McKernan et al J (120101) show that warm clouds from the surrounding halo occasionally impact galactic centers, 
potentially delivering lO'^-lO^M© of gas. 

Bondi-Hoyle-Lyttleton accretion causes gas to flow towards the black hole. At some point, the angular momentum 
of the gas will become important and the accreting gas must form a di sk. The final stage of accreti on is then governed 
by this accretion disk, which may be a geomet rically thin, radiative ( Shakura and Sunvaevlll973h or a geometrically 
thick, radiatively inefficient (e.g. ADAF) flow ( Naravan and Yil 1994 ). The details of the flow may be important for 
determining the spin of the black hole (see below) and the eff^ects of feedback from any nuclear activity (see 36.3.2l i. 

Galaxy-galaxy mergers can also lead to galaxies containing two (or potentially more) SMBHs, resulting in the 
potential for black hole mergers. The process of bringing two SMBHs together begins by dynamical friction against 
the background of dark matter (the same process wh ich is causing the blac k hole host galaxies to merge). The 



subsequent merging process was originally outlined bv lBegelman et al 
form a binary system at an initial separation of 



1980h and assumes that the two black holes 



flb 



G(Mi + M2) 

' 



(89) 



htt p : //lisa.nasa. gov/] 
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where Mi > M2 are the masses of the black holes and cr is the velocity dispersion of the host galaxy, such that 
the binary orbit contains a mass in stars and dark matter comparable to the sum of the black hole s masses. Around 
this b inary, the stellar distribution is expected to form a cusp with density profile p* oc r"^^'* dBahcall and Wold 



19761) . Initially, the binary hardens due to dynamical friction against the stellar background acting on each black 
hole individually. As the binary hardens this process becomes less effective as perturbations from distant stars tend 
to perturb the center of mass of the binary withou t changing its semi-major axis. However, once the binary becomes 



sufficiently hard, at a separation of dOuinlanll 19961) 



GM2 

Oh = (90) 
it can harden further by three-body interactions in which a passin g star is captured and then ejected at high velocity. 



The timescale, a /a, for hardening in this regime is dOuinlanll 19961) 



1 pa 



where the dimensionless hardening rate H ^ \5. If hardening continues long enough, the binary eventually becomes 
su fficiently hard that gravitational radiation dominates the evolution of the system which then coalesces on a timescale 



of dPetersll 19641) 



5c4fl4 

^'^^ ~ 256G3MiM2(Mi +M2)' ^'^^^ 

The ffy in the ointment of this neat picture is that the above estimates for the rate of hardening by stellar encounters 
assumes a fixed stellar background. In reality, as stars are ejected in three-body encounters the parts of phase-space 
containing stars that can be captured becomes depleted (so-called "loss cone depletion"). This inevitably slows the 
hardening process. The past ten years have seen numerous studies of this process and examination of various mech- 
anisms by which the loss cone may be refilled. For example, Yu ("2002!) finds that in triaxial potentials scattering 
of stars can efficiently refill the loss cone while iGo uld andRix (2000) suggest that the presence of a gaseous disk 
surrounding the black holes can help harden the binary. Numerical and analytical w orks have also indicated that 
the random walking of the binary center of mass may help mitigate loss cone depletion fouinl an and Hernquistlll99"7 : 



iMilosavljevic and Merritt 2001), while the effective refil ling of the loss cone by ejected stars returning on eccentric or 
bits (jlVIilosavlievic and Merridll2003 ; Sesana et al.ll2007 ) and interactions with stars bound to the binary ( Sesana et al 



I2OO8T may enhance the rate of hardening. While the details remain uncertain it seems that this basic process can lead 
to black holes merging in less than 10 Gyr. 

In addition to their mass, cosmological black holes are characterized by one other parameter, their angular mo- 
mentunQ The spin of a black hole can have a strong influence on the radiative efficiency and jet power of black holes 
(see ^6. 3. 2b and so the cosmological evolution of this quantity is important to understand. There are fundamentally 
two mechanisms which change the spin of a black hole: merging with another hole and accretion of material. 

The outcomes of binary black hole mergers have proven very difficult to simulate numerically. However, re- 
cent advances in n umerical techniques have allowed for successful simulation of the entire merging process (e.g 



Tichv and Marronettil2007l) and, theref ore, determination o f the final spin of the merger product. While the number of 



simulations carried out to date is small, iBovle et al. IfeOOS*) exploit symmetries of the problem to construct simple fit 



ting formula which accurately predict the spin of the final black hole as a function of the incoming black hole masses, 

spins and orbital properties^ 

As first considered by iBardeenI (Il970 l). material accreted from an accretion disk carries with it some angular 



momentum (approximately equal to its angular momentum at the latest stable circular orbit, before it began its plunge 
into the black hole). Defining a dimensionless spin parameter for a black hole through 

i - JclGMl (93) 



They will not possess any significant charge, as this would be quickly neutralized by accretion of oppositely charged material. 
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where J is the angular momentum of the hole of mass M,, such that < j < 1, then IShapird (l2005b defines a 
dimensionless spin-up function si f) by 

d/ M. 
df M._o 

where M, o is the rate of rest mass accretion. For a standard, relativistic, Keplerian thin disk accretion flow with no 
magnetic fields, the expected spin-up function due to accretion is given by (IShapirdl2005h 



= -Cisco - 2;£'isco, (95) 

where 

XiscoO) = , (96) 



nsco -\/'"?sco ~ 3risco + 2; yrisco 



and 



P 4co-2nsco + ;Vni^ 

£isco(;) = (97) 



'"isco ^r\^co ~ ^nsco + 2j Visco 



are the (dimensionless) specific angular momentum and specific energy of the innermost stable circular orbit (ISCO) 
of the black hole respectively. The radius of the ISCO orbit (in units of the gravitational radius GM./c^) is 

nsco = -V[3-Ai(;)][3+Ai(;) + 2A2(;)] 

+3+A2(;) (98) 



where 



AiU) = 1+[(1 -/)'''] 

x[(l+./)'/3 + (l-;y/3], (99) 

Alii) = ^3f+Ai(j)2. (100) 



Thornd (119741) found a small correction to this formula due to the fact that the hole preferentially swallows negative 
angular momentum photons resulting in a spin-up function for the standard thin disk that is positive for all j < 0.998, 
and therefore lets the black hole spin up to j ^ 0.998 in finite time as noted by Shapiro (2005). For a thick accretion 
flow the result is somewhat different (since the flow is no longer supp orted against gravity by its rotation alone as it has 
significant thermal and magnetic pressure). Be nson and Babul ( 20091) compute the spin-up function for ADAF models. 



Additionally, Benson and Babul (2009) compute how the m agnetic torques which allow black holes to drive jets (via 
the Blandford-Znajek and/or Blandford-Payne mechanisms; Blandford and Znaiekl l 1977tlBlandford and Pavnelll982h 



result in a braking torque on the black hole, spinning it down and resulting in an equilibrium spin of /' ^ 0.93 for a 
hole accreting from an ADAF. 

The relative importance of mergers and accretion for determining the spins of cosmological black holes depends 
upon the rate of galaxy mergers, the supply of gas to the black hole and additional factors such as the alignment of 
accretion disks and black hole spins and mergin g black hole spins and orb its. Man y of these factors are not too well 
understood. However, cosmological calculations dBerti and Volonteril2008t see also Volonteri et alJ2005 ) suggest that 



accretion dominates over mergers in terms of determining the spins of supermassive black holes, with the consequence 
that most such holes are predicted to be rapidly spinning. 

One additional consequence of black hole mergers is that the gravitational waves emitted during the final inspiral 
carry away linear momentum resulting in the black hole recoiling in the opposite direction with a velocity potentially 
large enough to unbind it from the galaxy (Fitchett 1983; Favata et al. 2004; Merritt et al. 2004; Blanchet et al. 2001). 
Studies of this effect within a cosmological framework suggest that such ejected black holes (wandering through 
intergalactic space) can make up 2-3% of the total mass density of supermassive black holes, while in individual 
cases they can account for up to half of the total black hole mass associated with a galaxy (.Libeskind et aL,2006.) . 
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1. 



6.3. Feedback 

As early as 1974 ( Larson 1974 : see also White and ReesI Il978l: bekeland Sill 1 19861) it was realized that star 
formation could not proceed with 100% efficiency in all dark matter halos. Evidence for this comes from multiple 
observed facts, but the two crucial ones are: 

The total mass density in stars is = (2.3 ± 0.34) x 10"^ (Cole et al.ll2001 ). much less than the total baryonic 
mass density of the Universe Qb = 0.0462 ± 0.0015 (iDunklev et al. 2009h . Therefore, only a small fraction of 
all baryons have been able to turn into stars. 

The distribution of galaxy luminosities (as described by the luminosity function) is very different from the 
distribution of dark matter halo masses. In particular there are many fewer faint galaxies relative to bright 
galaxies than there are low mass to high mass dark matter halos. If each halo conta ined baryons at the u niversal 
mix and turned all of them into stars we would expect these two ratios to be equal jBenson et al ] l2003h . 



This second point can be made even stronger: if all dark matter halos turned a constant fraction of their mass into stars 
we would still have too many faint relative to bright galaxies. There are a wealth of other observational constraints 
which indicate that the effic iency of galaxy formation must depend strongly upon halo mass — a useful summary of 
these constraints is given bv lBehroozi et al. I (l2010l) . 

Clearly what is needed is some process which preferentially suppresses the formation of stars in lower mass dark 
matter halos. The usual suspect for this process is energy/momentum input from supernovae explosions, perhaps 
augmented by stellar winds. 

A similar problem occurs in the most massive dark matter halos. Although cooling is relatively inefficient in such 
halos they can nevertheless cool significant mass of gas over cosmic time. Unchecked, this leads to the formation of 
galaxies significantly more luminous than any that are observed. The energetic requirements of this problem suggest 
that AGN are a possible solution. 



6.3.1. Supernovae/Stellar Winds 

Energy input from SNe due to a single stellar population are shown in Fig. [lOl If we assume that a fraction e of 
this energy is coupled into an outflow of ISM gas which leaves the galaxy at the escape velocity (so that we do not 
waste energy by giving the gas a velocity at infinity) this results in an outflow rate Mout given by 



1 . . 



- e 



Jo 



)£'sNe+winds(f - t')dt' 
M*(f)£'sN+winds, 



(101) 



where in the last step we have approximated energy input from SNe and winds as occurring instantaneously after star 
formation. This implies an outflow rate of 



Mo 



MJt) 



2e£sN+winds 

y2 ■ 



(102) 



2 

esc 



This has the required features: given some initial mass of cold gas, a fraction 1/(1 + /3) where p = 2e£'sN+winds/V, 
win be turned into stars, resulting in much lower star formation efficiencies in low mass galaxies. 

Figure[TO]shows the cumulati ve ener gy inpuj^into the interstellar medium as a function of time from a IMq burst 
of star formation with a Chabrieri (l2003h initial mass function. 

Similar arguments can be made using the momentum of supernovae explosio ns to drive the win d, which may 
be more relevant if energy is efficiently radiated from expan ding SNe bubble^ dMurrav et al 
momentum injection from supernovae is given by (IMurrav et al...2005i) 



2005). The rate of 



PsN ^ 2 X 10 
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kg m s 



(103) 



^*We assume that a Type II supernovae releases 10^' ergs of usable energy. We express the energy input in units of 10"ergs so that it can be 
interpreted as the equivalent number of supernovae. 
^'The momentum, of course, cannot be radiated away. 
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Figure 10: Cumulative energy input into the interstellar medium, expressed as the number of equivalent supemovae, per unit mass of stars formed 
as a function of time for a Chabrier initial mass function. The two panels show results for two different metallicities as indicated in the panels. The 
long dashed line indicates the contribution from stellar winds, the dotted line the contribution from Type la supemovae and the short dashed line 
the contribution from Type II supemovae. The solid line shows the sum total of all contributions. 



Since the momentum in any wind driven by this momentum input must satisfy 
velocity at infinity we have that 



A^fnufVco, where Voo is the wind 



Mo 



Vesc (Murray et al. 



(104) 

20051) . This gives a different scaling with the galaxy escape 



where we have assumed that 
velocity compared to energy driven winds. 

The details of how supemovae feedback actually operates remain somewhat unclear, and is complicated by the ne- 
cessity to understand how an expanding supemovae remnant interacts with a complex, multiphase interstellar medium. 
Our discussion below will mostly follow Efstathiou (2000; see also McKee and Ostrik er 1977). 

The basic picture considered involves a two phase interstellar medium consisting of hot, diffuse material (which 
fills most of the volume) and cold, dense clouds (which contain most of the mass). Evaporation of cold clouds by 
supemovae moves gas from the cold phase into the hot phase. Outflows, and suppression of star formation, from the 
hot phase will result if the hot phase: 

1 . fills most of the volume; 

2. is sufficiently low density that radiative cooling times are long, and; 

3. its temperature exceeds the virial temperature of galaxy (so that the gas may climb out of the potential). 

The effectiveness of feedback in such a model can be e xamined by consi dering the mass of cold clouds swept up 
and evaporated by an expanding supemovae remnant, M^y. lEfstathiou (l2000l) finds 



Me, 



(105) 



where £51 is the energy released by a supemovae in units o f 10^' ergs, Wh-2 is the numb er density of hydrogen in units 



of 10 cm , and the evaporation parameter as defined bv lMcKee and Ostriken(ll977h is 



2 = 



r 



47rfliA^ci0* 



(106) 



where y - Vh/ch ~ 2.5 relates the blast wave velocity, v-^, to the isothermal sound speed, Ch, ai is a typical cloud 
radius, A^ci is the number of clouds per unit volume, ipi^ relates the effective conductivity, /Ceff , to the classical thermal 
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conductivity, k, via /Ceff = (p^K and « 95 pc^ is the value of the evaporation parameter in the Solar neighborhood. 
(Units are chosen so that all variables have values of order unity in a typical galaxy.) 

Since Mgy is significantly larger than the mass of stars formed per supernovae then we may expect that if a 
significant fraction of this evaporated gas can escape as a wind then star formation may be efficiently suppressed. 

Gas in the hot phase will escape the galaxy if its specific enthalpy 

ivH^^ (107) 
2 2 p 

exceeds its specific gravitational binding energy. This therefore requires that we determine the thermodynamic prop- 
erties of the hot phase. Efstathiou (2000) does this by assuming that supernovae remnants have a porosity of unity (i.e. 
the total volume of remnants equals the volume of the galaxy, thereby creating an overlapping network of remnants). 
Given a model for the expanding supernovae remnant this allows the thermodynamic properties of the hot phase to be 
det ermined. 

Efstathioul (l2000l) shows that such a model of feedback can indeed produce self-regulated star formation in a 



galaxy disk that is quiescently forming stars (bursts are not required) and that both negative and positive (due to 
pressure-enhanced star formation) feedback can work simultaneously. 

6.3.2. AGN 

Very simple energetic arguments suggest that the energ y output involved in building a black hole in a galaxy can 



have significant effects on the formation of the galaxy itself dBenson et al .l2003h . To order of magnitude, let us assume 



that each galaxy contains a supermassive black hole with a mass equal to 0. 1 % of its stellar mass. The energy released 
by the formation of the black hole per unit mass of stars formed is then 

fi. = 10"^— ^c^ (108) 
1 - e 

where e is the efficiency of conversion of rest mass into energy output (eiffier in radiation or mechanical outflow). 
For typical values of e ~ 0.1 this implies & - 1.8 x 10^° ergs M^', which is an order of magnitude greater than that 
released by supernova explosions and stellar winds and is comparable to the energy released by cooling gas in forming 
some of the most massive galaxies. 

Feedback from AGN has the potential to directly link the properties of supermassive black holes and their host 
galaxies. As such, it may be a natural explanation for the observed correlations between supermassive black hole 
mass and galaxy mass or velocity disp ersion. The first study of how AGN feedback can lead to such correlation^ 



mass and galaxy mass or velocity disp ersion. Ihe tirst study or how AUJN reedback can lead to such correlationa i 
was described bvlsiik and ReesI (1199 8'). They showed that, on quite general grounds, coupling of the energy released 



by the formation of the supermassive black hole to the surrounding forming galaxy would lead to a relationship close 
to the observed M,-cr relation. The generality of these arguments imply that this result may be reasonably independent 
of the details of any more specific feedback model. 

The physical mechanism through which AGN feedback operates remains somewhat unclear AGN likely radia- 
tively heat cooling gas in ffie atmosphere surrounding galaxies, reducing the rate at which that gas can cool. Ad- 
ditionally, radiatively driven winds originating from the broad line region surrounding the black hole may result in 
mechanical feedback on the galaxy itself. Finally, in low accretion states, AGN may drive highly collimated and 
powerful jets which can reach out well into the surrounding halo. Plausibly all of ffiese mechanisms could be active 
during galaxy formatio n. 



Ciot ti et al.l (120091) examined the roles of radiative and mechanical feedback modes. They find that radiative 
feedback alone, while able to forestall the "cooling catastrophe" (see ^8.3l l. is unable to sufficiently limit the growth of 
black holes, resulting in them being too massive for a given galaxy (by a factor of around four). Mechanical feedback 
is found Jo be successful in limiting the growth of black holes. This is in a greem ent with other numerical simulations, 
such as lSpringel et al.i2005al:lMatteo et al.T2005l:ISiiacki et al.i2007i 120091: Ifiooth and Schavel2009i which incorporate 
energy injection from quasars (typically triggered by major mergers) and show that AGN activity can effectively 
expel gas from a galaxy and can establish the observed correlations between supermassive black hole and galaxy 



'Once again, assuming tliat any causative process is needed at all Jjahnke and Maccioll2010l) . 
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properties. However. ICiotti et al. I (l2009h find that mechanical feedback of high efficiency (as is required to match the 
scaling relations) is too effective, depleting galaxies of gas to a degree greater than is observed, while low efficiency 
mechanical feedback allows too much star formation to occur at late times, resulting in blue cores in elliptical galaxies 
that are inconsistent with observations. Precisely how efficiently jets from an AGN can couple their energy to the 
surrounding hot gas remains a topic of intense s tudy. The general p icture that is emerging is that the jets inflate bubbles 
or cavities in the hot atmosphere. For example. lOmma et al.l (120040 use numerical simulations t o show that jets inflate 



cavitie s in the hot gas that can excite g modes in the cluster atmosphere and heat the gas, while iRovchowdhurv et al 



(120041) explore a model in which the bubbles are buoyant and heat the atmosphere by doing pdV work as they rise. 

Recen t years have seen numerous implementations of detailed models of galaxy formation incorporating AGN 
feedba ck ( Scannapieco and OhI 2004; Croton et al. 2006 : Bower et al. 2006t Cattaneo et al. 2006 : Somerville et al 



2008a) . These have shown that AGN feedback can explain the long standing issue of the exponential break in 



the galaxy luminos ity function (without AGN feedback such models tend to produce far too many bright galaxies; 



Benson et al.l 120031) and helps produce a bimodal distribution of galaxy colors (as AGN feedback makes massive 
galaxies "red and dead"). These models typically a ssume a "radio mode" fe edback in which the AGN drives jets 



out of the galaxy while in a radiatively quiet state. iBenson and Babull (l2009l) give an expression for the jet power 
produced by a spinning black hole accreting from a geometrically thick accretion flowF^ relevant to this situation. As 
a result this type of AGN feedback is effective only when there is a quasi-hydrostatic atmosphere of hot gas in a dark 
matter halo for it to couple to. Therefore, AGN feedback only works in the virial shocking regime, resulting in it 
being important above a critical halo mass of order lO'^M© (with a weak redshift dependence). 



6.4. Chemical Enrichment 

The first generation of stars (known as Population III) and galaxies must have formed from primordial gas which is 
(almost) metal free. Stellar nucleosynthesis and subsequent pollution of the ISM and IGM (through stellar winds and 
supernovae explosions) with heavy elements has a significant impact on later generations of galaxies. In particular, 
the presence of heavy elements significantly alters the rate at which gas can cool (see ^3.2b and leads to the formation 
of dust which both attenuates optical and UV light from galaxies and re-emits that light at longer wavelengths. To 
accurately model the properties of galaxies therefore requires a treatment of chemical enrichment. 

The fraction of material returned to the ISM by a stellar population as a function of time is given by 



R{t) 



[M - M,{M\Z)]4>{M) — 

JM(i;Z) 



M 



(109) 



where <p(M) is the initial mass function normalized to unit stellar mass and M^(M) is the remnant mass of a star of 
initial mass M. Here, M(t) is the mass of a star with lifetime t. Figure[TT|shows the fraction of material recycled to the 
interstellar medium from a single stellar population as a function of time, f, since the birth of that population. Typical 
Pop. II initial mass functions (e.g. Salpeter 1955; Chabrier 2003) lead to around 40% of mass being recycled after 
10 Gyr, dropping to around 30% at 1 Gyr. A Pop. Ill initial mass function recycles much more mass (since it consists 
of much more massive stars), recycling 70% of the mass within a few tens of millions of years. 
Similarly, the yield of element ; is given by 



r'^" dMo 
V,(0= M,(Mo;Z)0(Mo) 

Jm(/;Z) 



Mo 



(110) 



where Mi(Mo;Z) is the mass of metals produced by stars of initial mass Mq. For a specified IMF we can compute 
yi(t, Z) for all times and elements of interest. Figure[T2]shows examples of the total metal yield for Pop. II and Pop. Ill 
initial mass functions. Typical total metal yields for Pop. II are around 0.04 increasing to 0.1 for the Pop. Ill initial 
mass function shown here. Stellar data are taken from lPortinari et aL ( 1998 ) for low and intermediate mass stars and 
lMarigoi ( i200ll) for high mass stars. 

In many cases, the process of chemical enrichment has been simplified by adopting the instantaneous recycling 
approximation, in which mass and metals are assumed to be returned to the interstellar medium instantaneously after 



'The jet power from a hole accreting from a thin accretion disk is much lower, due to the much weaker magnetic field in such a disk. 
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Figure 1 1 : The fraction of mass from a single stellar population born at time, ( = 0, recycled to the interstellar medium after time /. Solid and 
dotted lines show results for a Chabrier initial mass function with z ero and Solar m etallicity respectively, while the dashed line shows results for 
zero metallicity (Pop. Ill) stars with a lognormal IMF (case A from iTumlinsonl200(i) . 




Figure 12: The total metal yield from a single stellar population born at time, ( = 0, after time Solid and dotted hnes show results for a Chabrier 
initial mass function with ze ro and Solar met allicity respectively, while the dashed line shows results for zero metalUcity (Pop. Ill) stars with a 
lognormal IMF (case A from lTumlinsonl2006l) . 
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the formation of a population of stars (essentially assuming that stellar evolution happens instantaneously). Typically, 
the recycled fractions and yields are evaluated at some time, f, of order the age of the Universe. This approximation is 
reasonable at late times, when the ages of typical stellar populations in galaxies is several Gyr. For example. Fig. [12] 
shows that the metal yield is approximately constant after about 0.1 Gyr, although the recycled fraction (see Fig.fTTTi 
does not converge so rapidly. The advantage of the instantaneous recycling approximation is that it greatly simplifies 
the equations governing chemical enrichment such that we have 



JM(lia;Z) 



[M - M,iM;Z)](f>iM)^ s R 
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with a similar result for yields. The normal convolution integral to determine the rate of return of mass to the inter- 
stellar medium. 



then simplifies to 



MR(f)= R(t - t')M^(t')dt' 
Jo 

MR(f) = RM^it). 



(112) 
(113) 



This removes the need for knowledge of the entire star formation history of a galaxy to determine the rate of recycling 
(which can greatly reduce computational demand in numerical models of galaxy formation). 

While the instantaneous recycling approximation has its advantages, it is, of course, an approximation and one 
which breaks down in regimes of rapid star formation and at high redshifts where stars are much less than 10 Gyr old 
(since the Universe itse lf is very young). Detail ed model s of chemical e volution are well established wit hin the mono- 
lithic collapse scenario ( Matteucci and Greggial986: Matteucci and Gibs oal995; Timmes et al.. 1995; Franqois et al 



2004; Pipino and Matteucci 2004 20061 : Romano et al. 2005 ). but have only received limited study within the phys 
ically motivated context of hierarchical galax y formation. N umerous hydrodynamical s imulatio n codes now con- 
tain implementatio ns of chemical enrichment ( Scannapieco et al.. 2005 : Kobayashi et al.li2007; Pontzen et al. 2008 : 
Gnedin et al. I l2009h . S emi-analytic models of ^^ l ^^V formation are beginningjj to imple ment detailed chemical en- 
richment models also (Nagashima et al."2005a'b: ' Pipino et alj|2008l: lArrigoni et aLl l2009^. Within the context of hi- 
erarchical structure formation models such chemical enrichment modeling ha s been shown to be in a greement with 
measurements of the intrac luster medium metallicities of individual elements ( Nagashima et alj[2005ah . to reproduce 
the observed trend of metallicity and [ff/Fe] with stellar mass in elliptical galaxies (Arrigoni et al. 2 0091) and to repro- 
duce the measured distribution of metallicities in damped Lyman-a systems (iPontzen et al.ii2008i) . 



6.5. Stellar Populations 

Observational studies of galaxies make use of radiation emitted by them (or, sometimes, the lack of radiation due 
to absorption by galaxies) to infer their physical properties. As such, it is often crucial to be able to compute the 
emergent spectrum from each galaxy in a theoretical model. This is two-step process, involving first computing the 
spectrum of light emitted by all of the stars (and possibly the AGN) in the galaxy and, secondly, computing how this 
light is reprocessed by the gas and dust in the galaxy and along the Une of sight from the galaxy to the observer. 

6.5.1. Steller Population Synthesis 

In the absence of absorbing dust or gas, the luminosity as a function of frequency for a galaxy, the spectral energy 
distribution (SED), is simply a sum over the SEDs of its constituent stars. We can imagine a galaxy as consisting of 



Implementation of chemical enrichment models of this type is, in many ways, significantly easier in N-body simulations than in semi-analytic 
models. N-body simulations that include star formation typically spawn new star pailicles during each timestep to represent the stars formed at 
that time. Each particle can therefore be trivially tagged with a formation time and metallicity. On subsequent timesteps the rate of mass and metal 
return from that particle is easily computed and applied to surrounding gas particles for example. Semi-analytic models on the other hand have 
traditionally not recorded the full star formation history of each galaxy (due to computational resource requirements) as is needed to compute the 
chemical enrichment at each timestep. This can, of course, be done, as we are now seeing. 
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numerous "simple stellar populations" — ^populations of stars of the same age, metallicity and initial mass function. 
The luminosity, Ly^^''"'''\ of a galaxy at frequency v is then simply given by the following convolution integral 



r (galaxy) 



= I df' I dZ'M*(f',Z') 
Jo Jo 

x4^^P'(f-f',Z',0) 
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where M*(f, Z) is the rate of formation of stars of metallicity Z in the galaxy at time t and Lf^^\t, Z, (p) is the luminosity 
of a simple stellar population of age f, metallicity Z and with initial mass function <p{M). Given a model which predicts 
the rate of star formation in a galaxy as a function of time and metallicity (which is the galaxy formation aspect of this 
calculation), the problem is reduced to determining suitable (the stellar astrophysics aspect of the calculation). 

This, in turn depends upon the spectra of individual stars, Lv'^'*(f, Z), and the stellar initial mass function, (f>(M), such 
that 

Lf^^\t,Z,(f>)= 4>{M')Lf'''Ht,Z)m' , 
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where Mmin and M^ax are the minimum and maximum masses for stars respectively. 

r(SSP) 



Several libraries exist which provide L\, '(t, Z, (p) for several different ages, metallicities and initi al mass functions 



teruzual and Charlotll2003l: lMarastonll2005l: lOotter et alj|2007t IConrov et al.ll2008l; [Lee et al.ll2009l) . These are con- 
structed using a combination of theoretical stellar evolution models, observations of stars of known age and metallicity 
and theoretical models of stellar atmospheres where no good observations exist. While e normous progress ha s been 
made in understanding the spectra of stars significant uncertainties remain. For example. IConrov et al ] (2008), who 
performed a study of how uncertainties in such models propagate into constraints derived from galaxy observations, 
find that current models do not fully characterize the metallicity dependence of the thermally pulsating asymptotic gi- 
ant branch phase and that uncertainties in the slope of the initial mass functiorf^lead to an uncertainty in the evolution 
of the K-band magnitude of a stellar population of around 0.4 magnitudes per unit redshift (which leads to significant 
uncertainty when trying to interpret or predict the evolution of galaxy populations). 

6.5.2. Dust Absorption and Re-emission 

The presence of dust in galaxies has been known for a long time, and the effects of this dust on the observed 
properties of galaxies have been extensively studied. Simply put, dust absorbs light emitted by stars (and AGN), 
particularly at short wavelengths, is heated by this light and therefore re-emits it at longer wavelengths (typically in the 
infrared and sub-mm). The presence of dust in galaxies can therefore significantly affect their observed luminosities 
at optical and UV wavelengths. Models of galaxy formation must therefore take into account the effects of dust 
before comparing their predictions to observational data. The simplest such approach, adopted by many early models, 
is a "obscuring screen" or "slab" geometry, in which a plane of dust is placed in front of the model galaxy, given 
an extinction curve measured from the Milky Way (for example) and normalized to have an optical depth at optical 
wavelengths based on properties of the model galaxy (gas content and metallicity). This provides a simple estimate 
of the amount of extinction, but clearly does not reflect the true geometry of the dust (which is distributed throughout 
the galaxy) or the fact that stars may preferentially form in dense, dust regions. 

,Cole et al., (.2000,) introduced a model for dust extinction in galaxies which significantly improved upon earlier 
"slab" models. In Cole et al. the mass of dust is assumed to be proportional to the mass and metallicity of the 

ISM and to be mixed homogeneously with the ISM (possibly with a different scale height from the stars) and to have 
properties consistent with the extinction law observed in the Milky Way. To compute the extinction of any galaxy, a 
random inclin ation angle is select ed and the extinction computed using the results of radiative transfer calculations 
carried out bv lFerrara et al. I (ll999l) . 

Beyond these relatively simple models of dust extinction it has recently become possible to employ much more 
reaUstic ray tracing techniques to compute the effects of dust on galaxy spectra. For example, the Grasil software 



''See Kroupa (2001) for a recent discussion of the difficulties associated with determining the local initial mass function. 'Kroupa' l'200l'), along 
with Chary ( .2008) and .van Dokkum 12008,) . also discuss evidence for the non-universality of the initial mass function, although no definitive 
evidence for such yet exists. 
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Figure 13: Left-hand panel: An example of the spectral energy distributions from two face-on spiral galaxies as computed by Grasil iGranato et alj 
12000*) . The SED shows direct emission from stars at short wavelengths, with reprocessed light emitted fr om dus t at long wavelengths. In the 
intermediate region lines associated with PA Hs are visib le. Reproduced, with permission, from Granato et al] )2000l) . Right-hand panel: An image 
of a post-merger galaxy from a simulation bv lCox et alj flOOSi^ as co mputed by Sunrise iJonssor>i2006i) . A complex dust morphology, including a 
strong dust lane, is clearly visible. Reproduced, with permission, from lCox et alj i2008l) . 



dSilva et al.lll998l) is designed to compute the radiative transfer of star light through an ideaUzed galactic geometry 
consisting of a disk and a bulge each of which may contain both diffuse and clumpy gas and dust. Grasil takes an 
input galactic SED together with physical parameters of a galaxy (size of each component, metallicity and mass of 
gas present) and computes the resulting SED including absorption, scattering and remission from the dust, taking into 
account a realistic distribution of grains and polycyclic aromatic hydrocarbons (PAHs) the temperature distribution 
of which are computed self-consistently, and assuming that stars are born in the dense molecular clouds and escape 
from these on some timescale (leading to enhanced absorption in the UV which is produced primarily by young stars) 
(Fig.fTSll. This code is therefore ideal for (semi-)analytic studies of galaxies and has been employed to examine the ex- 
pected properties of submilhrneter and infrared galax ies in hierarchical cosmologies ( Granato et al J200(i Baughetal 
12004'. '2005; 'Lac ev et alJEoO S': 'Swinbank et al."2008') 



Similarly, the Sunrise code of Jonsson (2006) and Radishe bv lChakrabarti and Whitney (l2009h solve essentially 
the same problem (radiative transfer through a dusty medium) but work for arbitrary geometry using Monte Carlo, 
polychromatic algorithms and so are particularly well suited to hydrodynamical simulations of galaxy formation, 
allowing realistic images of si mulated galaxies to be made at any wavelength. Sunrise has been applied to studie s 
of luminous infrared galaxies (lYounger et al. I l2009h . quantitative morp hology of merger remn ants (ILotz et alJl2008l) . 
physical models to infer star formation rates from molecular indicators ( Naravanan et al. ''2010') and tests of our ability 
to recover physical parameters of galaxies from their broadband SEDs ( Wuyts et al. 2009). 

Since these ray tracing methods are computationally expensive they cannot yet be applied to large samples of 
model galaxies. The refore, attempts have b een made to construct simple r algorithms whic h capture most of their 



results. For example. iGonzalez-Perez et alJ (12008) extended the model of iCole et alJ (l200(]|) bv assuming tha t some 
fraction, /doud, of the dust is in the form of dense molecular clouds where the stars form (see lBaugh et alJ2005 ). Stars 
were assumed to form in these clouds and to escape on a timescale of Tquies ( for quiescent s tar formation in disks) 



or Tburst (for star formation in bursts), which is a parameter of the dust model dGranato et alj|2000 ). Since massive. 



short-hved stars dominate the UV emission of a galaxy this enhances the extinction at short wavelengths, so these 
stars spend a significant fraction of their lifetime inside the clouds. 
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To compute emission from dust they assumed a far infrared opacity of 



for A < /Ibreak 
forA>A break, 



(116) 



where the opacity normahzatio n a.\. A\ — 30 nm i s chosen to be a:i = 140cmVg to reproduce the dust opacity model 
used in Grasil, as described in Silva et alj (Il998h. The dust grain model in Grasil is a slightly modified version of 
that proposed by Draine and Lea i 1984^). Both the braine and Lee (1198 4') and Grasil dust models have been adjusted 
to fit data on dust extincti on and emission in the local ISM (with much more extensive ISM dust emission data being 
used by jSilvaetaL 1998). The normalization is set at 30/zm because the dust opacity in the IPraine and Led d 19841) 
and Grasil models is well fit by a power-law longwards of that wavelength, but not shortwards. The dust luminosity 
is then assumed to be 

Lv = 47r4v)Bv(r)Mz,gas, (117) 

where By{T) - [2hv^/c^]/[exp(hv/kr) - 1] is the Planck blackbody spectrum and Mz,gas is the mass of metals in gas. 
The dust temperature, T, is chosen such that the bolometric dust luminosity equals the luminosity absorbed by dust. 

This extended dust model, including diffuse and molecular cloud dust components, provides a better match to the 
detailed results of Grasil while being orders of magnitude faster, although it does not capture details such as PAH 
feature^; 

Fon tanot et al. l (l2Q09bl) have explored similar models which aim to reproduce the results of Grasil using simple, 
analytic prescriptions. They found that by fitting the results from Grasil they were able to obtain a better match to the 
extinction in galaxies than previous, simplistic models of dust extinction had been able to attain. 



6.6. Absorption by the Intergalactic Medium 

In addition to having to pass through the internal gas and dust in each galaxy, light emitted from a galaxy must pass 
through the entire intervening intergalactic medium between it and ourselves before we can observe it. The intergalac- 
tic medium contains significant amounts of neutral hydrogen (even at relatively low redshifts, long after reionization) 
which is known to be clumped into clouds by virtue of observations of the Lyman-o' forest of absorption systems 
in quasar spectra. Light emitted at redshift Zem at some wavelength /Igm, shortwards of the Lyman-a wavelength of 
1216A, will eventually be redshifted into the Lyman-c line and so will be absorbed by any clouds of neutral hydrogen 
at redshift 1 + Zabs = (1 + Zem)[/lera/1216A]. Light emitted at even shorter wavelengths may be additionally absorbed 
by higher lines in the Lyman series, or by the Lyman continuum shortwards of 912A. 

Models of the effective optical depth due to this absorpt i on (b ased upon the observ ed distribution of Lyman- 
a forest absorber properties) have been described by Madau ('1995') and Meiksid (l20Q6l) . The resulting absorption 
of starlight shortwards of Lyman-or is so severe that galaxies appear essentially dark at wavelengths shorter than 
/lobs — 1216A(1 + Zem)- This forms the basis of the "dropout" techniques for ide ntifying high redshift galaxies and 
quasars on the basis of their broadband colors (IMadau et al ] [l996t iFan et al ]|200lh . 



7. Computational Techniques 

The process of galaxy formation involves nonlinear physics and a wide variety of physical processes. As such, it is 
impossible to treat in full detail using analytic techniques. There are two major approaches that have been developed 
to circumvent this problem. The first, numerical N-body simulation, attempts to directly and numerically solve the 
fully nonlinear equations governing the physical processes inherent to galaxy formation. The second, semi-analytic 
modeling, attempts to construct a coherent set of analytic approximations which describe these same physics. Each 
has its strengths and weaknesses, as will be discussed below. As a result of this and of our ignorance of how some 
key processes work, both approaches contain a little of each other (rather like yin yang). A different, more empirical 
approach, utilizing so-called "halo occupation distributions" has become widely used in the past ten years and will be 
discussed briefly below. 
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7.1. N-body/Hydro 

The most accurate computational method for solving the physics of galaxy formation is via direct simulation, in 
which the fundamental equations of gravitation, hydrodynamics and perhaps radiative cooling and transfer are solved 
far a large number of points (arranged either on a grid or following the trajectories of the fluid flow). I will not attempt 
to review the numerical methods utilized in this approac h in any detail here (recent treatments of this can be found 



Bertschinger 1998 : Agertz et al. 2007 : Rosswoj 20091) . but will instead merely highlight some of the approaches 



used. 

Collisionless dark matter is (relatively) simple to model in this way, since it responds only to the gravitational 
force. For the velocities and gravitational fields occurring during structure and galaxy formation non-relativistic 
Newtonian dynamics is more than adequate and so solving the evolution of some initial distribution of dark matter 
(usually a Gaussian random field of density perturbations consistent with the power spectrum of the CMB) reduces 
to summing large numbers of forces between pairs of particles. In practice, clever numerical techniques (such 
as part icle-me sh, tree algorithms etc.) are usually used to reduce this A^^ problem into something more manage- 
able dKravtsov et al.lll997l: iSpringell 2005 ). The largest pure dark matter simulations carried out to date contain 



around 10 billion particles (*Spring el et al. 2005b|). Dark m atter only simulations of this type (carried out primar- 



ily for the cold d ark matter scen ario, but see White et all tl984: Klvpin et al. 1993; Bode et al. 2001; D ave et al 



20011: IColm et allbooi lAhn and Shapiroll2005h have been highly successful in determining the large scale struc 



ture of the Universe, as embodied in the so-called "cosmic we b". As a result, th e spatial and vel ocity correla- 
tion properties of dark matter and dark matter halos ^avis et al. lll985l; IWhite et al.l,1 987a b; Efstathiou etai]|l988l; 

Kravtsov et al.l 2004; Reed et al 



EkeetalJ[l996l; Ijenkins et al.1ll998l; IPadilla and Ba ugh 2002; B ahcaU et all 120041: 

2009), together with the density profiles (Navarro et al. 1997; Bul lock et al.P 2001b': 'Navarro et al."2004; Me m t et aL 
2005 JPrada et al. 20( )i), angular momenta (Barnes and Efstathio ulfl987 'lEfstathiou et al. 1988; WaiTen et al .l|l992 : 



Cole and Lacev 1 996tlLemson and Kaufftnann 1999; Bullock et al. 2001a; van d en Bosch e t al. 2002; Bett et al. 2007; 
Gao et al.tt200 8) and internal structure (lMoor~et al.„1999i:iKlvpin et al.il999i:iKuhlen et al.l2008i:iSpringel et al.i2008l) 



of dark matter halos are known to very high accuracy. 

Of course, to study galaxy formation dark matter alone is insuflicient, and baryonic material must be added in 
to the mix. This makes the problem much more diflicult since, at the very least, pressure forces must be computed 
and the internal energy of the baryonic fluid tracked. Particle-based methods (most prominently smoothed particle 
hydrodynamics; IS pringel 2005) have been successful in this area, as have Eulerian grid methods jRicker et al. 2000l; 



Frvxell et al.ll200Qc .Plewa and MiilleriuOOlt IOuihsll2004l) . The addition of radiative cooling is relatively straightfor- 



ward (at least while the gas remains optically thin to its own radiation) by simply tabulating the rate at which gas of 
given density and temperature radiates energy. 

Going beyond this level of detail becomes extremely challenging. Numerous simulation codes are now able to 
include star formation and feedback from supernovae explosions, while some even attempt to follow the formation 
of supermassive black holes in galactic centers. It should be kept in mind though that for galaxy scale simulations 
the real physics of these processes is happening on scales well below the resolution of the simulation and so the 
treatment of the ph ysics is often at the "subgrid" l evel, which essentially means that it is put in by hand using a semi- 
analy t ic approac h (Thacke r and Cou chman 2000; Kav et al. 2002; Marri and White 2003; White 2004: Stinson et d.l 
2006UCox et al.''2006: Scannapieco et al. 2006a; Tasker and Brvan 2006; Stinson 2007: Vecchia and Schava^OOsT 
Oppenheimer and Dave 2008; Okamoto et al. 2008b; Booth and Schaye 2009; Ciotti 200^^ 

Beyond this, problems such as the inclusion of radiative transfer or magnetic fields complicate the problem 
further by introducing new sets of equations to be solved and the requirement to follow additional fields. For 
example, in radiative transfer one must follow photons (or photon packets) which have a position, direction of 
travel and wavelength, and determine the absorption of these photons as they traverse the baryonic material of the 
simulation, while simultaneously accounting for the re-emission of absorbed photons at other wavelengths. De- 
spite these complexities, progress has been made on these issues using a variet y of ingenious numeri cal te chniques 



^el et al. 1999; Ciardi et al. 2001; Gnedin and Abel 2001; Ricotti et al. 2002; Petk ova and S pringel 2009; ' Li et al 
2008: Aubert and Tevssier.,2008!: .Colhns et al.. .2009; .Dolag and Stasvszva,2008; Finlator et al..,2009; Laursen et al 



20091) . Providing Moore's Law continues to hold true, the ability of simulations to provide ever more accurate pic- 



tures of galaxy formation should hold strong. A brief survey of current cosmological hydrodynamical codes and their 
functionalities is given in Table [1] 
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Table 1: A survey of physical processes included in several of the major hydrodynamical codes. The primary reference is indicated next to the name 
of the code. Where implementations of major physical processes are described elsewhere the reference is given next to the entry in the relevant 
row. 



Feature 


Gadget- 3' 


Gasoline^ 


HART^ 


Enzo(Zeus)'* 


Flash^ 


Gravity 


Tree 


Tree 


AMR'' PM ' 


AM#1P]VC1 


Multi-grid 


Hydrodynamics 


SPH^ 


SPtP 


amrP 


amrP 


AM# 


— > Multiphase subgrid model^ 




X 


N/A 


N/A 


N/A 


Radiative Cooling 


/ 


/ 


/ 


/ 




— > Metal dependent 




X 






✓El 


— > Molecular chemistry 


/15 


X 






X 


Thermal Conduction 




X 


X 


X 


/ 


Star formation 




/20 






X 


SNe feedback 






jm 


J2i] 


X 


— > Chemical enrichment 










X 


Black hole formation 




X 


X 


X 


/23 


AGN feedback 




X 


X 


X 


X 


Radiative transfer 


OTVET24.25 


X 


0TVETf2l 


/26 


/27 


Magnetic fields 


/28 


X 


X 


/29 


/30 



Notes 

"GAlaxies with Dark matter and Gas intEracT" JSpi-ingell . 1200^ : 

^Wadslev et aTl <2004h ; 

^Hydrodynamic Adaptive Refinement Tree jKravtsov et alil2002l) : 

^O'Sheaetal. (2004); 

'http://flash.uchicago.edu (' Frvxell et allbOOOl) : 

''Adaptive Mesh Refinement; 

'Particle-mesh; 

**Smoothed Pailicle Hydrodynamics; 

'Applicable only to SPH codes — used correctly, AMR codes naturally resolve multiphase media; 
Scannanieco et al. 1 2006a); 
"Banerjee et al. (2006); 
' iScannapiecp et al 1 200^ ; 

Tassis et airfoosl) ; 
'iSmithetalJhoOgl); 
''Yoshida et al. (2^; 
"'Equilibrium only; 
''Turk (2009); 
'^■Jubelgas etaDfe004); 
''Scannapieco et al. (20o3); 
2'^Govemato et al. (2007); 
2'Tasker and Bryan (200^; 
^^Matteo et al .. t200a); 
2'Federrath et alj <2010l) ; 
^''Optically Thin Variable Eddington Tensor; 

2'Petkova and Springel (2009); 

^''Flux-limited diffusion approximation ( No rman et al J206^ ; see also lWise and Abell2008bl) ; 
2'Riikhorst et al. ( 2006); Peters et al. (201(J); 

^^Dolag and Stasyszyn (2008); 

"Collins etal. (2 009; see also .Wang and Abel2009l) ; 
^' lRobinsonetalJj2004l) ; 
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7.2. Semi-Analytic 

The technique of "semi-analytic modeling" or, perhaps, "phenomenological galaxy formation modeling" takes 
the approach of treating the various physical processes associated with galaxy formation using approximate, analytic 
techniques. As with N-body /hydro simulations, the degree of approximation varies considerably with the complexity 
of the physics being treated, ranging from precision-calibrated estimates of dark matter merger rates to empirically 
motivated scaling functions with large parameter uncertainty (e.g. in the case of star formation and feedback — just as 
in N-body/hydro simulations). 

The primary advantage of the semi-analytic approach is that it is computationally inexpensive compared to N- 
body/hydro simulations. This facilitates the construction of samples of galaxies orders of m agnitude larger than 
possible with N-body techniques and for the rapid exploration of parameter space (IHenriques et al. 2009) and model 
space (i.e. adding in new physics and assessing the effects). The primary disadvantage is that they involve a larger 
degree of approximation. The extent to which this actually matters has not yet been well assessed. Comparison studies 
of semi-analytic vs. N-body/hydro calculations have shown overall quite good agreement (at least on mass scales well 
above the resolut ion Umit of the simulation) but have been limited to either s implified physics (e.g. hydrodynamics 
and coohng o nly; Ben son et al .l2001bl:lYoshida et alJ2002HHellv et alJ2003bl) or to simulations of individual galaxies 
dStringeretal . 2010). 

Some of the earliest a ttemp t s to co nstru ct a self-consist e nt sem i-analytic model of galaxy for mation began with 
White an d Frenk" (mi), ' Cold d 199 lb and iLacev and SiUd (Il99lh, dr awing on earlier work by 'Rees and Ostrikeil 
d 19771) an d White and Re esI (Il978l)7since then numerous studies jKauf fmann et al. 1993; Baugh et al. 1999bl ll998 
Somerville and Primackll999l;ICole et allbOOOtlBenson et alJl2002al;lHatton et alj2003l:lMonaco et alj2007l) have ex 



tended and improved this original framework. A recent review of semi-analytic techniques is given bv lBaughl(l2006t) . 
Semi-analytic models have been used to investigate many aspects of galaxy formation including: 



37Lh 
ughl 



galaxy counts dKauffmann et alJ[l994 : Devriendt and Guiderdonill2000f) : 



galaxy clusteri ng dOiafe rio et alJI 19991: Kauffmann et aDll999allbl : iBaugh et allll999at Isenson et al.ir2000cllal : 
Wechsler etal] l2001: Bla izot et alJbOOd)" 



galaxy colors and metal licities dKauffmann and Chariot! 1998l:ISt?ringel et alJ200lt iLanzoni et al.l2005l:lNagashima et al, 
2005bl:lFontetaDl2008l) : 



sub-m m and infrared galaxies dGuiderdoni et alj 119981 : iGranato et al. I l2000t iBaugh et al.l 120051 : iLacev et al 
2008h : 



abundance and properties of Local Group galaxies ( Benson et alj|2002bt Somervill3l2002 : Li et al ] l20ll : 
the re ionization of the Universe ( Devriendt et aljl998 ; Benson et alj2001a ; Somerville and Liviol2003 ; Benson et al 



2006h : 



the heating of galactic disks dBenson et al ] |2004 : 

the properties of Lyman-break galaxies dGovernato et al ] |l998l:lBlaizot et alJl2003ll2004l:lLacev et"ai]l2010l) : 



supermassive black hole formation and AGN feedbac k (Kauf fmann and Haehnelt [20001: ICroton et alJ 1 120061: 
Bower et alj|2006t iMalbon et al1l2007l: ISomerville etalj,2008b: iFontanot et aL,2009'l : 



damped Lyman-a systems dMaller et al ] |200lll2003h : 

the X-ray properties of galaxy clusters dBower et al ] |200lll2008h : 

chemical enrichment of the ICM and IGM dLucia et al ] |2004lNagashima et alj|2005al) : 



the formation histories and mo rphological evolution of galaxies dKauffmannll996HLucia et alj2006l: l iFontanot et al 
2007tlSomerville et al.lE008al) . 
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The term "semi-analytic" model has become somewhat insufficient, since this name now encompasses such a 
diverse range of models that the name alone does not convey enough information. Semi-analytic models in the 
literature contain a widely disparate range of physical phenomena and implementations. Table |2] is an attempt to 
assess in detail the physics included in currently implemented semi-analytic modelQ 



73. Halo Occupation Distributions 

Given the complexity of galaxy formation it is sometimes desirable to take a more empirical approach to the prob- 
lem. This can be advantageous both to relate o bservations to the assumed underlying physical structure of the Uni- 
verse (e.g. IZheng et al.ll2007t [Tinker et alJl201Q) . and to make predictions based on extrapolations from curr ent data 
and w hich have a cosmological underpinning. "Halo occupat ion distributions", first described by|N evman and Scott 
( 1952 ) and used to study galaxy clustering (see, for example, Benson et al. 2000c; Peacock and Smithll2000l) in their 
simplest form specify the probability of finding galaxies of some prescribed type in a dark matter halo of mass 
M, P(N\M). Given this probability distribution, knowledge of the distribution of dark matter halo masses and their 
spatial distribution plus some assumptions about the locations of galaxies within halos, one can construct (statistically 
or within an N-body simulation) the distribution of galaxies. This allows, for example, the abundance and clustering 
properties of those galaxies t o be inferred. A detailed discussion of the use of halo occupation distributions in studying 
galaxy clustering is given by Coorav and ShethI ( 2002h . 

This approach has recently been extended by incorporating simple prescriptions relating, for example, star for- 
mation rates or quasar activity to halo mass and redshift. Several authors have demonstrated that such empirical 
models, while simple, can fit a wide variety of galaxy data and can be use d to gain insight into phenomena such as 
"downsizing" dYang et all l2003t IConrov and Wechslen , 120091: 1 ICrotonll2009i) . At a fundamental level, the success of 
simple, empirical models such as these suggests that, despite the complexity of galaxy formation physics, its outcome 
is relatively simple. 



8. Topics of Current Interest 

Having reviewed the major inputs to our current theory of galaxy formation and explored some of the tools 
employed to solve them, in this section we briefly explore five topics of particular interest in contemporary galaxy 
formation theory. 



8.1. The First Galaxies 

After the initial excitement of the Big Bang, inflation, nucleosynthesis and the recombination the Universe enters 
a protracted period of relative quiescence, known as the "Dark Ages", during which there are no luminous sources 
and structure grows primarily only in the dark sector. The end of the Dark Ages is thought to occur when the first 
star forms. The formation of these first stars is an unusual problem in astrophysics in that it represents a clean, well- 
defined problem that can be tackled numerically — the initial conditions are simply those of an expanding Universe 
containing small perturbations in the dark matter and baryonic components which is a s imple mi xture of hydroge n and 
helium. As such, this problem has been examined and explored in significant detail dAbel et alj|2002l : lYoshida et alj 
20061 : Gao et al. 2007 ). The consensus picture that has emerged from this work is that the first star will form due to 



gas coUapse and cooling (via the gas-phase formation of molecular hydrogen) in a halo of mass around 10 at a 
redshift of z ~ 50. The first star is likely to be very massive (current simulations cannot follow the evolution of the 
proto-star beyond the point at which it becomes optically thick in various lines of molecular hydrogen due to a lack 
of treatment of the necessary radiative transfer, but at this stage the collapsing core typically has an accretion rate that 
is expected to lead to the formation of a star with a mass of several tens of Solar masses) and, therefore, short lived. 
What happens after this first star forms rapidly becomes a much more complicated problem, involving the poorly 
understood physics of supernovae explosions and the numerically challenging problem of radiative transfer. 



^^One component missing from Table|2]is tlie specific implementation of stellar population synthesis used by each code. We have chosen to not 
include this because all the models listed treat this calculation in precisely the same way (see ^6.5. IK differing only in which compilation of stellar 
spectra they choose to employ. Since it is trivial to replace one compilation of stellar spectra with another we do not consider this a fundamental 
difference b etween the models. Neve rtheless, the choice of which stellar spectra to use can have important consequences for the predicted properties 
of galaxies iConrov and GunnLl2010l) and so should not be overlooked. 
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Table 2: A survey of physical processes included in major semi-analytic models of galaxy formation. In each case we indicate how this process 
is implemented and give references where relevant. In many cases a single model has implemented a given physical process at different levels of 
complexity/realism. In such cases, we list the most "advanced" implementation that the model is capable of. 



Feature 

Merger Trees 

— > Analytic 

-> N-body 
Halo Profiles 
Cooling Model 

— > Metal-dependent 
Star Formation 
Feedbacks 

-> SNe 

-> AGN 

— > Reionization 
Merging 

Substructure'^ 
Substructure-Substructure'' 
Environments 

— > Ram Pressure Stripping 

— > Tidal Stripping 

— > Harassment 
Disks 

Disk Stability 
— > Dynamical Friction^^ 
— > Thickness 
Sizes 

Adiabatic contraction 
Chemical Enrichment 
Dust 



Durham' 

Modified ePS'' 

Einasto'^ 

/ 
/ 

/ 

/'2 
JM 

N-body 1^ 

^/20fl0l 



/ 
/ 

/ [delayecPl 
Grasil^^ 



Munich 

ePS^ 
/ 

Isothermal 

/ 
/ 

/ 
/ 
X 

N-bod>El 

X 

/24 

X 
X 

/ 

X 
X 

X 

/ [instant^'] 
Screen^^ 



Model 

Santa-Cruz^ Morgana^ 

ePS PiNoccHio** 

/ X 

NEW NEW 



/ 
/ 

/ 
/ 
/ 

^21,22 

X 
/ 
X 

/26 
X 
X 

/ 

/ [delayed^"] 



Slab 



,34 



/ 
/ 

/ 

/ 

DlGH 

X 
X 

/ 

X 

/ 

X 
X 

/ 

/ [instant] 
GrasilE2135 



Galic? 

X 

/ 

Empirical" 

/ 
/ 

/ 

N-bod>El 
/ED 

/ 
X 

/ 

X 
X 



/ [delayed^'] 
SlalP 



Notes 



fCole et al.' 12000'); 
'Croton et al. (200'3); 
jSomerville et al. ( 20083); 
''Monaco et al 12007) ; 
'Hattoneta lJ<2003l); 

^Parkinson et al.l J20Q^; 

' Kaulfmann and White 0993h ; 
* Monaco et al. 120020; 
•^ Helly et al. 1,2003^ 



iBenson and Bowd <2010h ; 



A "dark mat t er" co re is included in calculations of disk sizes with an empirically selected dark matter fraction; 
^Bower et al.' ^200^); 
' Cattaneo et al. (20M); 
"Maccio et al. l20(m 
' Lanzoni et al. 1 2003); 

'How does the model track substructures within halos?; 

'Substructure orbits and merging times are determined from N-body simulations; 

**Dynamical Friction; substructure merging times are computed from analytic estimates of dynamical friction timescales; 

'Does the model allow merging between pairs of subhalos orbiting in the same host halo?; 
^"Using hierarchically nested substructures; 
^'Using random collisions of sub halos; 

Somerville and P rimackl il999l) ; 49 
2'*Fontetal. (2008); 
^''Briigeen and Lucia 1 200^; 
^'Lanzoni et al. (2005); 
2'' Somerville et al. 1 2008a); 

^ 'Poes the model incl ude dynamical friction forces exerted by a galaxy disk on orbiting satellites?; 
^ ^iBenson et"ai] )20Q4) ; 
^''Lucia et al.'l'2004); 
Arrigoni et al. 1 200^; 



3C 

^ipipino etal. l200t 
^^Silvaetal. (1998); 



Needless to say then, studying the formation of the first galaxies is more complicated still, although no more 
complicated than modeling the formation of later generations. In the next decade new facilities such as JWST and 
thirty meter class telescopes will begin to open a window on the earliest generations of galaxies and, as such, there is 
a need for theoretical understanding and predictions for this epoch of cosmic history. In particular, this represents a 
real opportunity to make testable predictions from galaxy formation theory which can be confronted with data in the 
near future. 

In considering the earliest generations of galaxies there are various physical processes which, while often neglected 
at lower redshifts, are of potentially crucial importance. For example, understanding the prop erties of P opulation III 
stars, their various feedback processes (b oth positiv e and negative) and initial mass func t ion (iGlover and AbeL2008 

3 

may De c rucial tor unde rstanding how this popul; 



Trenti and Stiavelli| 120091: lOhkubo et alJi 2009: Ma chida et alJ 120081: iLawlor et alJ 120081: lO'Shea and NormanlEoO? 
20081 : Harris et al. 2007; Stacv andBromm 2007) may be c rucial for unde rstanding how this populati on of galax- 
ies grows CTegmarket al..,1997- Ricotti et al.. 2001.: Iwise and Abel2007a.bl: lOreif et al.ll2008l: IWise and Abel.2008b; 



Wise et al. 2008; Bromm et al 



20091) . Additionally, connecting these galaxies t o later generations involves under- 



stand ing how they evolve through (and cause) the epoch of cosmic reionization (IBenson et al.ll2002at lAlvarez et al 



20081) . modeling non-local feed backs (i.e. interac tions between neighboring galaxies or the galaxy population as a 
while) and feedbacks in general jBusha et al ] l2009l) . 

Various interesting questions await definitive answers. For example, how long is the period of Population III star 
formation? This seemingly simple question depends on how we define Population III — theoretical reasoning suggests 
that we may expect a significant shift in the initial mass function once gas is raised above a critical metallicity of 
IO^^Zq (Sa ntoro and S hull 2006) and the ability of Population III stars to enrich their surroundings (and, therefore, 
future generations of star formation). Since Population III stars are likely to produce many more hydrogen ionizing 
photons per unit mass than their Population II counterparts this question may have important consequences for the 
reionization history of the Universe. 

Reionization itself remains an interesting problem. While it has been convincingly demonstrated theoretically 
( Chiu and Ostriker 2000t Ciardi et al. 2003 : Somerville et al. 2003 : Benson et al. 2006h that t he universe could be 
reionized by galaxies sufficiently early to match constraints on the optical depth to the CMB (iDunklev et al. I I2OO9I) 
these calculations are currently forced to make assumptions about the distribution of mass and ionized regions in the 
IGM and about the ionizing photon escape probability from high redshift galaxi es ( Dove et al.l2000t Ricotti and Shull 
2OO0I : iRazoumov and Sonmier-Larsen 200(51: Gnedin et al. 2008; Yaiima et al. 2009t Wise and Cen 20091) which are 
still poorly understood. 

In addition to observing high redshift galaxies directly, it may be possible to learn much about them from the 
study of their present day remnants . "Galactic arche ology" o f this sort has b een investigated quite extensively 
(IFreeman and Bland-Hawthorn 20021; IScannapieco et al.t,2006b: Brook et al ] |2007h . 

The conclusions are that Population III stars (usually taken as a convenient proxy for "first galaxies") that formed 
in progenitor galaxies of the Milky Way likely formed over a fairly broad range of redshifts, with the formation rate 
peaking at z ^ 10 but continuing until z ~ 4-5 (in chemically isolated halos which have not been contaminated by 
metals produced by star formation in earlier forming halos). As such, any remnants of these stars are probably spread 
over most of the nearby Galactic halo. The lack of detection of such metal-free stars therefore places constraints on 
their initial mass function — they must have been sufficiently massive to have evolved off' of the main sequence by 



the pr esent day. Currently, this suggests that such stars must have been more massive then about O.8M0 (iBrook et al 



20071). If we examine truly old stars, however, the situation is rather different. Simulations suggest that these formed 
in halos close to the high density peak that eventually formed the final system, and therefore are preferentially located 
in the central regions (typically the bulge) of the Milky Way. 



8.2. The Formation and Sizes of Galaxy Disks 

Galactic disks are, perhap s , the most prominent feature of galaxies in the local Universe. Additionally, recent 
observations dLaw et al. 1 120091: ISchreiber et alJ I I2OO9I) have demonstrated that massive disks are already in place as 
long ago as z = 2. Understanding the details of how disks form and grow is therefore of crucial importance for galaxy 
formation theory (ISilldl200lh . 

The basic picture has been in place for a long time (iFall and Efstathioul[l980l) . Collapsing systems of dark matter 
and gas are imparted some small amount of angular momentum as a result of tidal torques from inhomogeneities 
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in the surrounding matter distribution (iHovle [1943 iBarnes and Efstathioulll987n . While gas can radiate energy and 
thereby collapse to form a galaxy it cannot radiate angular momentum. As such, angular momentum is conservecQ 
during the collapse and is eventually responsible for halting that collapse when the proto-galaxy reaches a radius at 
which its angular momentum is sufficient to provide support against the gravitational potential in which it resides. 
Back of the envelope type estimates o f the resulting sizes of disks show that th is picture produces disks with sizes 
comparable to those that are observed dFall and Efstathioul[l980l: lMo et alj|1998l) . However, this basic picture makes 
several simplifying assumptions, some of them difficult to justify. A careful analysis of these assumptions is therefore 
warranted to understand the process of disk formation in greater detail. 

It has been convincingly established that the distributi on of angular momentu m of gas in halos prior to any cooling 
and galaxy formation is similar to that of the dark matter ( van den Bosch et alJ20 02. 2003) in so far as the distribution 
of spin par ameters is very similar (the directions of the angular momenta vectors of dark matter and gas are often 
misaligned: Ivan den Bosch et alj|2002h . The distribution of the angular momentum of g as within an individual halo 
is less well studied, but recent hydrodynamical simulations ( Sharma and Steinmetzll2005 ) have demonstrated that the 
distribution of the specific angular momentum, j, can be well described by a simple form: 



1 dM _ 1 
M d7 ~ f^T(aj) 



(118) 



where T is the gamma function, M is the total mass of gas, /'d - jtot/o: and jtot is the mean specific angular momentum 
of the gas. The parameter aj takes on different values in different halos, with a median value of 0.89. The fraction of 
mass with specific angular momentum less than j is then given by 



fi<j) = y\aj,j-j. 



(119) 



where y is the incomplete gamma function. Despite this knowledge of the initial distribution of angular momentum in 
the baryonic component, the question of how well angular momentum is conserved during the process o f cooling and 



collap se is less clear However, recent simulation studies have begun to shed some light on this issue. IZavala et al 



(l2008h describe high resolution SPH simulations in which the baryonic physics (particularly the star formation and 
feedback from supernovae) is varied such that one simulation results in the formation of a disk galaxy and another 
in the formation of a spheroidal galaxy. They show that, in the case of the disk galaxy the angular momentum of 
the baryonic material that will eventually form the galaxy tracks the behavior of the dark matter angular momentum: 
growing as predicted by linear theory at early times and then remaining constant. In the case of the spheroidal galaxy, 
the angular momentum still grows initially as expected from linear theory, but then 90% is rapidly lost as pre-galactic 
fragments undergo mergers which transfer their angular momentum to the dark matter In the disk galaxy case, strong 
feedback prevents gas from cooling into these small halos at early times, and so avoids angular momentum loss to the 
dark halo. 

In addition to macroscopic quantities such as the size of a disk, the picture described above can predict the radial 
density profiles of disks. This is interesting since observed disks typically follow an exponential radial profile (at least 
out to some radius, beyond which the stellar light profile typically truncates — gaseous disks typically extend out fur- 
ther beyond this stellar truncation radius). In a hierarchical Universe, in which dark matter halos and their gas content 
grow in mass and size with time, the angular momentum of material accreting onto a galaxy will, on average, increase 
with time also. As a result, disk formation is an "inside out" process, with early infall leading to the formation of the 
inner regions of the disk and later infall adding primarily to the outer regions. Knowing the distribution of specific 
angular momenta of the gas from which the disk formed it is possible to determine the distribution of radii at which 
that gas will settle and, consequently, the radial density profile. Of course, this assumes that the angular momenta of 
gas elements are unchanged after they arrive in the disk. In reality this may not be the case — viscosity in a gaseous 
disk may redistribute angular momentum and even in a purely stellar disk stars can scatter from perturbations such 
as spiral density waves and thereby change their angular momentum and may help to establish the exponential radial 



It is conserved if we ignore the possibility of torques between baryons and dark matter which could transfer angular momentum from the gas 
to the dark halo. We will review this assumption later 
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profile if the viscous timescale is of ord er the star formation timescale ( Silk and Normanl l 198ll:lLin and Prinridl 19871: 
Clarkdl 19891 1991uFirmani et alj|199 6). Current results suggest that exponential disks can plausibly arise given our 
present understanding of the initial distribution of angular momentum, s tar formation efficiency and supern ovae feed- 
back without the need for viscous redistribution of angular momentum ( Duttonll2009t Stringer et al. 20 id) , although 
the detail ed properties (e.g. scale length distributions) of the disks do not necessarily agree perfectly with what is 
observed (IDutton et al.luOOTh . Of course, comparison with measurements of stellar disks requkes determination of 
not just the mass density profile but the stellar density profile. Star formation will proceed most rapidly where the gas 
is densest (i.e. in the inner regions of the disk) and may be inhibited altogether at larger radii where the background 
radiation can prevent the formation of molecular hydrogen and where the disk may become stable to perturbations 
and so unable to form molecular cloud complexes (iMartin and Kennicuttll200ll:ISchavell2004 . 

The assumption that the angular momentum distribution of the disk remains constant after formation is not pre- 
cisely correct even in the absence of hydrodynamical effects such as viscosity. Stars in a disk can interact with per- 
turbati ons such as spiral density waves, exchanging energy and angular momentum with them. ISellwood and Binnev 
(l2002h show that interactions between stars and spiral waves at the corotation resonance lead to an exchange of energy 
and angular momentum that changes the radius of a star's orbit while keeping it circular A star can, in principle, in- 
teract with multiple spiral waves and thereby migrate over a sig nificant radial range. The question of how much radial 
mixing occurs due to this effect has been examined recently bv lRoskar et al. ( 2008a ) and Roskar et al. (2008b!)- They 
find that this is a significant process which can lead to the formation of a population of stars beyond the star formation 
truncation radius in a disk, and which also flattens metallicity gradients (particularly in older stellar populations). 

Disks are dynamically rather fragile systems and so not only must they form but they must survive across cos- 
mological spans of time. This requires that the disk be stable to the growth of large scale perturbations which could 
destroy it. Additionally, in a hierarchical universe, survival is a non-trivial feat as the disk must survive frequent 
merging with other galaxies and constant bombardment by dark matter substructures orbiting within its own dark 
matter halo. Significant work has gone into studying the effects of dark matter substruc t ures on galactic diskj^ 
dToth and Ostrikeilll992l: IVelazauez and White! 1 19991: Ipont et al.ll200ll: iBenson et al.1 12003: iKazantzidis et al.1 120081 



Read et al ■I l2008t IKazantzidis et al.ll2009l: IPurcell et alJ l2009'). The current consensus seems to be that the cold dark 



matter model is consistent with the presence of thin galactic disks, but only just so, and some fraction of disks have 
probably been significantly thickened or destroyed by interactions with dark matter substructures. 



8.3. The Overcooling Problem 



White and ReesI (119781) proposed that galaxies would form as gas cools inside of dark matter halos and demon- 



strated that this provided a reasonable estimate of the typical mass scale of galaxies. This simple picture has a long 



standing problem however The mass function of dark matter halos rises steeply (approximately as M Reed et all 
2007h at low masses. Since cooling is very efficient in these low mass halos we might expect the galaxy mass and/or 
luminosity function to show a si milar slope at the l ow mass/luminosity end. In fact, measured slopes are much shal- 
lower (typically around -1; e.g. ICole et al.ll200lb . Rectifying this discrepancy is usually achieved by postulating 
some form of feedback, typically from supernovae, which ca n inhibit star f ormation in these low mass systems (e.g. 
by driving a wind out of the galaxy). However, as shown by Benson et al. (l2003h this causes another problem — too 
much gas is now left over to accrete into massive halos (group and cluster mass objects) at late time wherein it cools 
and forms over-massive galaxies (much more massive than any galaxy observed). 

This "overcooling" problem is not easy to solve, for the simple reason that the energy scales involved are much 
larger than for lower mass systems (the characteristic potential well depth of a dark matter halos scales as M^^^ for 
halos of the same mean density). Several possible solutions have b een proposed howe ver, ranging from t hermal 



conduction (Benson et al.ll2003t iDolag et al.ll2004t iPope et al.ll2005h . massive outflows (IBenson et al.l 120031) . mul- 



tiphase cooling (in which gas cools but is lock ed up into clouds which are inefficiently transferred to the galaxy; 
Mailer and BuUockll2004t iKaufmann et al ]l2009l) and feedback from active galactic nuclei. 

The AGN feedback scenario has gained considerable favor in the past few years for a variety of reasons. Firstly, 
observations have indicated that all galaxies seem to contain a central supermassive black hole with a mass that scales 



Primarily because this is potentially a means by which to rule our the cold dark matter scenario; if the predicted abundance of substructure 
was sufficient to destroy galactic disks, the model could be ruled out. 
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roughly in proportion to the mass of the galaxy dMagorrian et alJI 1998c iGebhardt et alj|2000at iFerrarese and Merritt 



2OOOI: IGebhardt et alj I l2000bl: iTremaine et al.ll2002l) . The formation of these black holes must have involved the lib 



eration of large amounts of energy and material sinks deep into the potential well of the black hole in some form of 
accretion flow (IBenson et al.ll2003h . If this energy can be successfully utilized to counteract the overly rapid cooling 
of gas onto the galaxy it would provide a natural source of sufficient energy available in every galaxy. Additionally, 
some feedback loop connection between supermassive black hole and galaxy formation of this type is attractive as it 
provides a means to explain the correlation between galaxy and black hole properties. 

Feedback from AGN can be divided into two categories: quasar mode and radio mode. The quasar mode is trig- 
gered when large amounts of gas are funneled into the circumnuclear disk from much larger scales by, for example, a 
merger between two galaxies in which torques act to transfer angular momentum from the gas, causing it to flow in- 
wards. This is likely the dominant mechanism for black hole mass growth and results in a high (relative to Eddington) 
accretion rate — most likely through a thin accretion disk — and significant optical luminosity. In contrast, radio mode 
feedback occurs when the black hole is accreting at a more modest rate due to Bondi accretion from a diffuse hot atmo- 
sphe re of gas and is in an optically dim but radio loud phase. In this phase, the black hole is thought to dr ive powerful 
jets jMeiejl 19991 12OOII: IVilhers et all 120051: iNemmen et aDl2007t iKomissarov et al.ll2007t iBenson and Babul 20091) 
which can reach to large distances and have been seen to have a significant impact on their surroundings (iBirzan et al 



20041) . The mechanism via which energy from the jets is efficie ntly coupled to the hot atmosphere of gas remains 



poorly understood: a combination of observational evidence (e.g. lOwen et al.ll2000) and theoretical insights suggest 
that jets inflate bubbles or cavities in the hot atmosphere (see ^6.3.21) . Despite these uncertainties, semi-analytic treat- 
ments which simply assume an efficient coupling have demonstrated that this can effectively shut down cooling in 
massive halos, resulting in a reduction in the masses of the largest galaxies and good agreement with luminosity func 



tions and the bimodal distribution of galaxy colors dCroton et 



2st galaxies ana good agreement witn luminosity tunc - 
all booel: iBower et alj|2006l: ISomerville et al.ll2008bh . 



Observational evidence in support of AGN feedback is beginning to emerge ( Schawinski et alj 2007h and"seems to 
favor a radio mode scenario (■Schawinskiii2009i) . 



8.4. Local Group Dwarf Satellites 

As with any theory, galaxy formation theory is often best tested by exploring its predictions in extremes of physi- 
cal conditions. For example, while our current theory may do well at explaining the properties of galaxies of average 
luminosity, it may perform less well in explaining the properties of the lowest mass galaxies that can form. Unfor- 
tunately, low mass galaxies are also low luminosity and so are difficult to observe. Fortunately, our own privileged 
position within the Milky Way and the Local Group means that there exists a population of faint, dwarf galaxies much 
closer by than we might otherwise expect. Study of the population of Local Group galaxies therefore has the potential 
to teach us much about galaxy formation at the lowest mass scales (and therefore, in a hierarchical Universe, at the 
earliest times). We must be cautious, however, of over-interpreting conclusions drawn from the Local Group which 
is, after all, just one small patch of the Universe. 

Until the advent of the Sloan Digita l Sky Survey around 1 1 satellite galaxies were known to exist within the 



Until tne aavent ot tne iiloan Digita l aKy aurvey arouna ii satellite galaxies were Known to exist witnin tne 
virial radius of the Milky Way's halco ( Mated 1998 ). The deep and uniform photometry of the SDSS has allowed 



many new such galax ies to be found by searching for overdensities in the stellar distribution in both posit ion and color 
JWillman et alJ2005h . The current roster of dwarf satellites within the Milky Way's halo amounts to 24 (ToUerud et alJ 
20081) . However, the SDSS covers only around one fifth of the sky and detection algorithms are imperfect. Re- 



cent estimates of the total number of dwarf satellites in the Milky Way's halo alone are of the order a few hundred 
jToUerud et alj 20081). The faintest galaxy currently known. Segue 1, has a luminosity of only 340 times that of the 
Sun (Beloku rov eTal]l2 007). 

These dwarf galaxies are interesting from a theoretical perspective for several reasons. F irstly, their low mass 



makes them highly sensitive to feedback effects, both internal (from supernovae for example; Font et al. 2010l) and 
external ( photoheating of the IGM by the entire population of galaxies can, in principle, inhibit the formation of these 
satellites; iGnedin l2000t lOkamoto et alj I l2008ah . Furthermore, they may have gained at least some of their mass as 
a result of molecular hydrogen cooling, and have very low metallicities which may provide insight into the process 



^'We are following the definition of virial radius given bv lBenson et all i2002al) here. Under this definition, 9 Local Group satellites are within 
the virial radius of M3 1 , while a further 1 9 lie outside of the virial radii of the halos of both the Milky Way and M3 1 . 
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of chemica l enrichment. Beyo nd these galaxy-related aspects, the dwarf satellites appear to be highly dark matter 
dominated dStrigari et al ]l2008h . making them excellent systems in which to explore the properties of dark matter and, 
therefore, to place interesting constraints on the behavior of dark matter on small scales. 

Consequently, these galaxies have attracted significant theoretical attention over the past few years. In particular, 
when simulations demonstrated that dark matter substructure could survive within a host halo for significant periods 
of time, it became clear that the Milky Way's halo should contain thousands of small dark matter halos that could 
potentially have formed a dwarf galaxy within them dMoore et al .ll999HKuhlen etalJ2008l: ISpringel et al .I2OO8I) . This 
was in stark contrast to the 10 or so satellites known at that time. This "missing satellites" problem appears to have 
been largely solved. Theoretical models have demonstrated that a combination of supe rnovae feedback, suppression of 



- J r-^r^ 1 

galaxy formation by a photoionizing background and tidal plus ram pressure stripping (IBullock et al J2000tlSomerville 
' 2002l:lBenson et alj|2002at iBusha et alj|2009l: iMaccio et alj|2009l: IWadepuhl and Springel"2QlQ' 'Font et al."2010') 



reduce the number of satelhtes significantly (leaving a population of truly dark subhalos; Kormendv and FreemanI 



20041) . while the discovery of new satellites (and accounting for ones as yet undetected) has raised the observational 



target closer to the theoretical predictions. 

Beyond this first order question though, there are many other interesting aspects to the Local Group satellites. For 
example, while they all c ontain populations of old stars, many of them show evidence of periodic star formation at 



much more recent times ( Tolstoy et al. 20091) . What triggers this star formation? Given that feedback can strongly 



affect these systems, it is interesting to ask if we can gain insight into the feedback mechanism by studying these 
galaxies (e.g. the ir metallicities ar e highly sensitive to the strength of feedback so may place useful constraints on any 
feedback model; Font et al. 2010). Finally, some of these galaxies must have been torn apart by the tidal field of the 
Milky Way (the Sagittarius dwarf is currently suffering t his fate) — the dwarf satellite population may therefore have 
contributed to the build-up of the Milky Way's stellar halo ( Cooper et alj2009l) . and there may be remnant tidal streams 



which could provide information about th e hierarchical growth of the Milky Way's halo ( Bullock and Johnstonll2005 



Robertson et aLi2005i : .Font et al...2006a.bc iJohnston et al.i,2008.) . 



8.5. The Origins of the Hubble Sequence 

The variety of morphologies of galaxies is, perhaps, the most obvious observed characteristic of galaxies. Tra- 
ditionally, morphology has been measured "by eye" by a trained observer who classifies each galaxy into a different 
morphological class based upon (amongst other things) the pr ominence of an y central bulge, how concentra ted the 
light distribution is and the presence or otherwise of dust lanes. Hubble ( 1936t see also^d e Vaucouleursll 19591) placed 
galaxies into a morphological class ification scheme using such an approach, and this basic morphological classifi- 



cation has persisted to the present (Ide Vaucouleurs et al. 1991). Applying this type of morphological classification 



to tod ay's large datasets is difficult, but has been achieved by utilizing "crowdsourcing" techniques (iLintott et al 



20081) . A key observational goal is to assess how the Hubble sequence evolves over time, as this should place strong 



constraints on theoretical models. This is observationally challenging, but progress is being made (e.g. Kriek et al] 
E009i) . 

Theoretically, morphological evolution is inevitable in a hierarchical universe in which galaxies interact with each 
other Early work demonstr ated that major mergers bet ween galaxies could transform disks into spheroids leading 
to morphological evolution (IBarnes and Hernquistlll992) . but it is clear that th e picture is more c omplicated, with 
secular evolution playing an important role in transforming disks into spheroids (iKormendv and Ken nicutt 2004) and 
major mergers of very gas rich systems can lead to t he reformation of a disk after the merger is over (lBarnesll2()02l: 
Springel and Hernquisli2005l; [Robertson et al.ll2006al) . 



A fundamental difliculty in assessing the ability of any given theoretical model to explain the morphological 
properties of galaxies is that the definition of morphology itself is very complicated, and somewhat nebulously defined. 
This problem is beginning to be circumvent ed, both by the ability of numerical simulations of galaxy formation to 
produce realistic "mock images" of galaxies dJonsso nl2006h which can, in principle, be classified by ey e just as a real 



galaxy, and by the use of more quantitative morphological measures suc h as bulge-d isk decompositio n dSchade et al 



19961: iRatnatun^a et al .Ul999tlSimard et al.U2002UBen son et allbOOTh. Gini-M20 (iLotz et alH2004l) . concentration- 
asymmetry-smoothness (Watanabe et al., 1985; Abraha m et al r il996HBellet Ill l200i) ' and so on. 



Understanding the build-up and evolution of the Hubble sequence of morphologies is, nevertheless, an impor- 
tant task for galaxy formation theory. The morphological structure of a galaxy clearly tells us something about its 
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formation history and therefore captures information which its stellar population, for example, does not. Several 

J Si — « — 1 ! ll 1 I 

Studies have explored this aspect of galaxy formation theory d Baugh et al. 1996a b; Kauffmann 1996; Govern ato et al 



1999: Firman i and Avila-ReesdEool iLucia et aP 120061: iLucia and Blaizot 2007: Governato et al. 2007; Pa rrv et aL 
20091: Benson an d"Devereuxll2009l) ^ with the general conclusion that hierarchical cosmologies can plausibly give rise 
to the observed mix of morphological types, although with significant uncertainties remaining in both the modeling 
of morphological transformation and in the comparison with observed morphologies. The consensus opinion is that 
massive elliptical galaxies form most of their stars in smaller progenitor galaxies and only assemble them into a single 
galaxy significantly later, while spiral galaxies are dominated by in situ star formation. Secular evolution of disks 
into spheroids also seems to be a crucial ingredient, particularly in the production of lower mass spheroidal systems 
dParrv et al.ll2009HKormendv et al.ll2009l) . 



9. Future Directions 

The fundamental goal of galaxy formation studies is to comprehend how the laws of nature turned a Gaussian 
random distribution of density fluctuations laid down by inflation into a complex population of galaxies seen at the 
present day. At this time, this author does not see any convincing evidence that any new physics is needed to explain 
the phenomena of galaxieQ The problem is more one of complexity: can we tease out the underlying mechanisms 
that drive different aspects of galaxy formation and evolution. The key here then is "understanding". One can easily 
comprehend how a 1/r- force works and can, by extrapolation, understand how this force applies to the billions of 
particles of dark matter in an N-body simulation. However, it is not directly obvious (at least not to this author) 
how a 1/r^ force leads to the formation of complex filamentary structures and collapsed virialized objects. Instead, 
we have developed simplified analytic models (e.g. the Zel'dovich approximation, spherical top-hat collapse models 
etc.) which explain these phenomena in terms more accessible to the human intellect. It seems that this is what we 
must strive for in galaxy formation theory — a set of analytic models that we can comprehend and which allow us to 
understand the physics and a complementary set of precision numerical tools to allow us to determine the quantitative 
outcomes of that physics (in order to make precision tests of our understanding). 

The division of galaxy formation models into N-body/hydro and semi-analytic is rather idealized. In reality there 
is significant overlap between the two — many semi-analytic models make use of dark matter halo merger trees drawn 
from N-body simulations while many hydrodynamics simulations include recipes for star formation and feedback 
which are semi-analytic in nature. It seems most likely that these two techniques will continue to develop and may, in 
fact, grow less distinguishable, incorporating aspects of each other into each other (once again, yin yang). 

Galaxy formation benefits from a wealth of observational data, often to the degree that it is an observationally 
lead field in which theory plays the role of trying to explain observed phenomena. This deluge of data is unlikely to 
cease any time soon — we can expect more and higher quality data on local galaxies and also the arrival of usefully 
sized datasets of galaxies at the highest redshifts. Galaxy formation theory should continue to attempt to develop our 
comprehension of these observed phenomena but should also strive to move into the regime of ma king true predictions 
for as-yet-unobse rved regions of parameter space (e.g. the high redshift, z >6, Universe; e.g. Finlator et al. 2010l; 



Lacev et al. hoiol) . Only in this way can we grow our confidence that we have truly understood the physics of galaxy 



formation. 



10. Summary 

Our theory of galaxy formation is gradually becoming more and more complete, but it is clear that large gaps in 
our understanding remain. This is, perhaps, not suprising — galaxy formation incorporates a wide array of physical 
processes, many of which we can currently observe the consequences of only indirectly. Most of the physics of 
galaxy formation is inherently nonlinear, making it difficult to obtain accurate solutions. Finally, what we observe 
from galaxies is usually several steps removed from the underlying physical properties (mass, density etc.) that we 
would ideally like to know about. Despite all of these difliculties, rapid progress is being made. The next decade 



By "new physics" here I mean modifications to established physical laws, new forces or fields etc. Of course, dark matter and dark energy 
probably require "new physics" of one type or another, but I will leave those as a problem for cosmology. . . 
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should see this trend continuing via a combination of ever better observational datasets (both low and high redshift) 
and the continued development of novel theoretical tools. 
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